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rends  in  augmented  system  dynamics  (eigenvalues)  were  shown  to  agree  with 
previous  results.  Moving-base  simulation  results  obtained  for  a  highly- 
banked  flight  condition  were(usedf  tojal so) establish  a  new  multi-axis  pilot 
model.  This  model  is  considered  suitable  for  simultaneous  control  of  the 
dynamically  interacting  longitudinal  and  lateral -directional  axes  of  the 
vehicle/display  system.  The  system  augmentation,  based  on  this  model,  again 
was  predicted  to  significantly  improve  performance.  The  trends  in  system 
eigenvalues  for  various  augmentation  levels  showed  good  agreement  with  the 
above  longitudinal  results.  However,  the  lateral-directional  axis  was  found 
to  be  much  more  important,  and  significant  modification  of  the  dynamics  were 
predicted  to  lead  to  system  improvements. 
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SECTION  I 


INTRODUCTION  AND  SUMMARY 

OBJECTIVES  OF  THE  RESEARCH 

With  the  advent  of  active  control  technology  (ACT)  and  control 
configured  vehicles,  the  use  of  specifications  on  the  rigid  body  modes 
of  an  aircraft  for  flight  control  design  becomes  unacceptable.  Reducing 
structural  weight  by  employing  active  flutter  suppression  or  stability 
augmentation,  and  enhancing  maneuverability  by  direct-lift  and  side- 
force  control  surfaces  are  some  of  the  topics  of  current  flight  vehicle 
research.  Such  vehicles,  however,  intrinsically  involve  dynamics  of 
higher  order  than  the  rigid  body  modes,  and  these  dynamics  have  been 
shown  to  significantly  alter  piloted  vehicle  performance  as  well  as 
subjective  pilot  ratings  of  the  systems.  The  inclusion  of  such  additional 
dynamics  in  the  determination  of  the  acceptable  open-loop  systems 
characteristics  so  drastically  increases  the  dimensions  of  the  problem 
as  to  render  this  approach,  impractical.  Furthermore,  this  is  attacking 
the  problem  "upside  down  and  backwards." 

The  ultimate  design  objective  is  to  maximize  the  performance  of  the 
man- vehicle  system.  The  most  logical  approach  would  seem  to  be  to  develop 
and  apply  systems  analysis  methods  based  on  closed-loop,  task-oriented 
design  techniques  -  representing  the  man  in  the  system  appropriately  as 
a  dynamical  decision  and  control  element.  Attacking  the  problem  this  way 
naturally  leads  to  "optimal  control  configured  vehicles"  and  forces  one  to 
ask  the  right  questions  along  the  way.  For  example,  we  naturally  address 
the  questions  of  what  physical  parameters  the  pilot  is  monitoring  and 
attempting  to  regulate  in  specific  piloting  tasks.  Only  after  answering 


such  fundamental  questions  can  one  make  an  intelligent  design  decision 
with  regard  to  factors  that  would  affect  these  rel(a^ejit  physical  para¬ 
meters. 

The  objective  of  this  research  activity  then  was  to  explore  the 
utility  of  multi-variable  control  techniques  for  piloted-flight-vehicle 
control  synthesis,  with  the  above  factors  in  mind,  by  application  to  a 
specific  piloting  task. 

The  methodology,  developed  previously,  employs  an  optimal  multi- 
variable  control  model  (OCM)  of  the  pilot,  specifically  includes  the 
pilot's  mission  objectives  for  control  synthesis,  and  simultaneously 
determines  the  pilot's  analytical  representation  (or  model)  and  augmenta¬ 
tion  control  law.  As  a  result,  the  method  is  intended  to  be  well  suited 
for  high-order  dynamical  systems,  inherently  reflect  the  cooperative 
structure  of  the  pilot  and  flight  control  system,  and  at  the  same  time 
include  the  pilot's  ability  to  adapt  to  the  controlled  system  dynamics. 

The  original  development  of  this  method  is  reported  in  Reference  1,  and 

the  approach  could  be  considered  an  optimal  control  version  of  a  control 

[21 

design  approach  utilizing  Anderson,  "paper  pilot".  J 
SUMMARY  OF  RESULTS 

The  piloting  task  addressed  throughout  the  study  was  air-to-air 

tracking.  The  initial  phase  involved  the  analysis  and  vehicle  augmentation 

for  single  axis  pitch  tracking,  the  system  dynamics  including  active  heads- 

up  display  dynamics  as  well  as  those  of  the  vehicle  rigid  body.  Fixed-base 

m 

simulation  data  obtained  elsewhere1-  J  were  used  to  determine  the  pilot 
model  parameters,  specifically,  the  objective  function  and  those  determining 
the  observation  and  neuromotor  noise  characteristics.  It  was  significant 
that  a  unique  objective  function  was  found  in  the  study  cited  for  the  task 
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in  question,  over  a  variety  of  vehicle  dynamics.  Such  an  assumption 
would  be  necessary  for  use  in  augmentation  synthesis. 

By  parametrically  varying  the  level  of  optimal -control  augmentation, 
not  only  rms  tracking-error  performance  improvements  were  predicted,  but 
trends  in  desirable  system  dynamics  were  obtained.  These  trends  were 
depicted  as  open-loop  as  well  as  piloted  closed-loop  system  root  loci 
and  the  agreement  was  noted  between  these  loci  and  the  asymptotic  behavior 
of  linear  quadratic-optimally-controlled  systems.  Significant  in  this 
regard  is  the  appearance  of  "pseudo  zeros"  or  transmission  zeros  in  the 
results  and  their  importance  in  interpreting  these  results.  Specifically, 
p  of  the  augmented  system  roots  tend  to  move  toward  the  transmission 
zeros,  and  the  n-p  remaining  roots  ultimately  approach  a  Bu^terworth-type 
pattern,  where  n  is  the  order  of  the  system  and  p  is  the  number  of 
transmission  zeros. 

The  dynamics  of  the  displays  were  found  to  affect  the  "optimal" 
vehicle  characteristics.  Agreement  was  found  between  the  results  for  the 
less  "active"  of  the  two  sights  and  previous  results  obtained  with  a  fixed 
sight  (and  no  display  dynamics).  However,  with  the  more  dynamic  of  the 
two  displays,  a  different  set  of  desired  vehicle  roots  were  obtained. 
Clearly  then  the  complete  system's  dynamics  must  be  considered,  not  just 
the  rigid-body  eigenvalues.  Furthermore,  the  numerator  dynamics  (which 
include  the  display  dynamics)  were  found  to  affect  the  desired  system 
characteristics  (e.g.,  handling  qualities). 

With  the  augmentation  determined,  the  pilot  model  was  then  used  to 
predict  rms  performance  (e.g.,  tracking  errors)  as  well  as  pilot  and 
vehicle  describing  function  frequency  response.  This  analysis  indicated 
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the  potential  for  significant  reduction  in  tracking  error  and  pilot 
workload,  with  an  antirioated  improvement  in  pilot  rating. 

The  second  phase  of  the  study  involved  a  similar  analysis  of  air-to- 
air  tracking  in  a  highly-banked  turning  flight  condition.  As  such,  it 
involved  pilot/vehicle  analysis  of  a  flight  condition  in  which  the 
longitudinal  and  lateral -directional  axis  were  not  independent,  but 
involved  unsymmetrical  cross  coupling.  Consequently,  a  significant 
amount  of  pilot-modeling  activity  was  necessary  as  no  multi-axis  pilot 
models  have  been  developed  and  documented. 

To  accomplish  this  modeling  activity,  results  from  a  large-amplitude 
motion  simulation  (the  Air  Force  LAMARS  facility  at  Wright-Patterson 
AFB,  Ohio)  were  compared  to  predicted  rms  performance.  It  was  determined 
that  in  addition  to  tracking  error  and  displayed  target  lead  angle,  the 
weighted  outputs  in  the  single-axis  pilot  objective  function,  the  pilot 
was  also  attempting  to  minimize  lateral  acceleration  and  relative  bank  angle 
between  the  attacker  and  target.  Hence,  the  pilot’s  vector  of  observations 
as  well  as  the  objective  function  weights  were  modified  appropriately. 

The  use  of  this  pilot-modeling  approach  for  tracking  a  target  not 
undergoing  random  maneuvers  was  another  unique  result  of  this  phase. 

Without  a  random  forcing  function  (or  process  noise)  the  pilot's  observation 
and  motor  noise  are  the  only  source  of  errors.  Since  this  was  actually 
the  situation  simulated,  the  modeling  proceeded  appropriately.  The 
solution  involved  selecting  the  variances  of  the  additive  motor  noises 
yielding  the  best  match  of  statistical  results  (i.e.,  rms  performance), 
and  a  result  of  this  modeling  task  is  the  noise-to-signal  ratio  for  the 
motor  noises. 

This  multi-axis  pilot  model  having  been  established,  the  multi- 
variable  augmentation  synthesis  followed  directly.  Again,  a  family  of 
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full-state  feedback  control  laws  were  synthesized  with  the  pilot-optimal 
control  formulation,  and  the  open-  and  closed-loop  root  loci  used  to 
depict  the  trends  in  system  dynamics  (in  terms  of  the  eigenvalues).  The 
longitudinal  results  were  found  to  agree  with  the  single-axis  results 
using  the  same  display  dynamics,  and  were  explainable  in  terms  of  the 
asymptotic  behavior  of  optimally  controlled  systems. 

The  more  complex  lateral  directional  axis  was  also  found  to  yield 
results  consistent  with  the  asymptotic  properties.  The  dutch-roll  mode 
was  augmented  to  yield  significantly  higher  frequency  and  damping, 
ultimately  approaching  a  Butterworth  configuration,  while  the  spiral  and 
roll  subsidence  were  coupled  to  form  a  complex,  stable  roll/spiral  mode. 

It  is  significant,  in  closing,  that  the  pilot/vehicle  analysis 
and  augmentation  results  all  point  to  the  high  relative  importance 
and  complexity  of  the  lateral  directional  axis.  The  parameters  associated 
with  this  axis  required  much  higher  attentional  allocation,  as  determined 
by  the  pilot  model,  and  the  dynamics  were  highly  modified  via  augmentation. 
Yet  this  axis  has  received  far  less  attention  in  the  literature  and 
even  less  experimental  activity  is  apparant.  Much  further  work  is 
therefore  called  for  in  our  opinion. 

ORGANIZATION  OF  THE  REPORT 

The  following  section  (Section  II)  includes  the  complete  analysis 
and  results  from  the  first  phase  of  the  study,  that  involved  with 
single-axis  pitch  tracing.  The  multi-axis  pilot/vehicle  analysis  and 
augmentation  is  reported  in  Section  III. 

A  significant  amount  of  the  effort  of  the  project  was  devoted  to 
modeling  the  vehicle,  display,  and  target  kinematics  and  dynamics,  these 
results  are  documented  in  Appendices  A  through  D.  The  matrices  defining 
the  linear  state  variable  model  of  the  multi-axis  tracking  problem  is  given 
in  Appendix  E. 
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Finally,  a  management  summary,  citing  publications,  students  involved, 
and  interactions  with  Air  Force  Laboratory  personnel  appears  in  Appendix  F. 


) 
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SECTION  II 


PITCH  TRACKING  ANALYSIS  AND  AUGMENTATION 

The  problem  to  be  addressed  is  the  analysis  and  augmentations  of  the 
plant  dynamics  in  an  (air-to-air  pitch)  tracking  task.  The  higher-order 
system  dynamics  actually  include  linearized  flight  vehicle  rigid-body  and 
actuator  dynamics  as  well  as  the  display  dynamics  of  two  lead-computing  sights. 
In  each  of  the  cases  considered,  the  vehicle  flight  condition,  and  therefore 
rigid-body  dynamics,  remain  unchanged.  However,  the  sight  dynamics  vary 
considerably  between  cases  with  different  sight  types  and  different  tracking 
ranges. 

The  vehicle  and  one  of  the  sights  was  identical  to  one  of  the  cases 
addressed  experimentally,  as  well  as  analytically  (via  the  OCM),  by  Harve^^ In 
his  investigation,  Harvey  used  fixed-base  simulation  data  to  infer  the  pilot's 
objective  function  in  this  task.  A  unique  objective  function  was  found  to 
yield  extremely  good  correlation  (in  terms  of  rms  statistics)  between  ana¬ 
lytical  and  experimental  results  over  a  wide  range  of  system  dynamics.  It  is 
therefore  assumed  in  our  work  here,  that  this  same  objective  function  is  invar¬ 
iant  with  the  system  dynamics,  and  is  used  in  all  cases  considered. 

The  objectives  of  this  part  of  the  investigation  then  include  not  only 
the  estimation  of  the  performance  improvement  attainable  with  augmentation,  but 
also,  and  perhaps  more  important,  the  determination  of  the  most  desirable  plant 
or  controlled  element  dynamics.  By  parametrically  varying  the  level  of  augmen¬ 
tation,  we  are  able  not  only  to  establish  trends  in  the  desirable  plant  char¬ 
acteristics  with  varying  high-order  plant  characteristics.  Finally,  the  results 
so  obtained  are  compared  to  those  of  Hollis^lho  used  simple  pitch-rate  and 
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plunge-acceleration  feedback  along  with  Anderson's  paper  pilot  to  determine 
optimum  vehicle  dynamics  for  pitch  tracking  with  a  fixed  sight  (with  no  dis¬ 
play  dynamics). 

THE  SYSTEM  MODELS 

The  vehicle  dynamics  are  the  linearized  short-period  approximation  for 
a  typical  fighter  aircraft  in  level  flight  at  15,000  ft.  altitude,  Mach  = 

0.9  (V  =  952  ft/sec),  and  the  specifics  are  presented  in  Table  1.  Also,  to 
aid  in  the 

TABLE  1 

VEHICLE  MATH  MODEL 

“  =  Za“  +  q  +  Z<55E 
0  =  q 

q  =  (Ma  +  M-Za)a  +  (Mq  +  M^q  +  (M^  +  M-Z^Sj, 

'65stick 

Za  =  -983.  (ft/sec2)  Ma  =  -10.4  (rad/sec2) 

Zfi  =  -90.5  (ft/sec2)  M^  =  -0.344  (rad/sec) 

Ta  =  .05  sec  Mq  =  -0.738  (rad/sec) 

Cs p  =  0.32;  u>sp  =  3.35  (rad/sec)  Mg  =  -37.1  (rad/sec2) 

interpretation  of  results,  we  will  consider  augmenting  the  simple,  second- 
order  plant  (e(s)/6st(s)  =  K/s2). 

The  schematics  of  the  displays  are  shown  in  Fig.  1.  The  tracking  task 
of  course  involves  the  desire  to  minimize  the  tracking  error  e,  where  for  the 
simple  K/S2  plant 


Display  Schematics 


e  =  e  -  8c 

For  the  pitch  tracking  task 

e  =  A  -  8 

where  8,  the  relative  1 ine-of-sight  angle  is  governed  by  the  relations 

8  =  q  -  q  (yj  -  e  +  a) 

yT  =  "  l  aT 


Note  for  this  display,  the  fixed  reference  mark  represents  the  weapon  line  and 
the  sight  (reticle  and  pipper)  is  dynamic  (moves  in  display). 

For  the  high-order  pitch-tracking  task,  the  displayed  variable,  or  lead 
angle  A,  includes  the  lead  required  for  velocities  and  acceleration,  as  well 
as  a  ballistic  "jump"  correction.  The  governing  linearized  equation  of  an 
ideal  sight  is  (from  Appendix  B) 

yT 

a  =  Tf(e  -  e)  -  20“  aj  ’  ~ET  a  0) 


This  relation  implies  that  the  1 ine-of-sight  rate  B  and  the  target's  normal 
acceleration  a^,  are  available  for  lead-angle  calculation,  and  therefore  rep¬ 
resents  an  idealized  sight.  A  second  lead-angle  equation  representing  an  actual 
typical  sight  is  given  as  (from  Apoendix  B) 


i  1  Tf 

=  .  1  »_/_S L_L 
Tf  A  '  2D 


.038  q)  a+0 


(2) 


This  relation  is  obtained  from  the  above  by  assuming  the  1 ine-of-sight  rate  6, 
may  be  approximated  by  the  lead-angle  rate  A,  and  the  target  acceleration  is 
approximately  equal  to  the  attacker's,  or  a^  -  la.  Note  that  since  the  actual 
sight  (Eqn.  2  )  depends  entirely  on  the  attacker's  variables  (a,e,  etc),  it  is 
much  more  sensitive  to  pilot  stick  input  than  the  ideal  sight.  Also  note  the 
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sight  dynamics  are  clearly  a  function  of  tracking  range.  Two  tracking  ranges 
were  considered,  0  =  1000  ft.  and  D  =  3000  ft,  associated  with  two  projectile 


times  of  flight,  Ty  =  0.33  sec.  and  Ty  =  1.30  sec,  respectively.  The  process 
driving  noise  for  the  pitch  tracking  case  is  the  target  normal  acceleration 
generated  by  filtered  white  noise  where 
ay  =  -.33ay  +  n 


n  =  -0.33n  +  w(t) 


and  the  white  noise  intensity  is  selected  to  yield  chosen  levels  of  target 

2 

accelerations.  The  driving  noise  for  the  K/s  plant  is  simply  the  commanded 
variable  ec,  generated  by  filtered  white  noise  where 

e  +  ao  +  be„  =  w(t) 
c  c  c 

For  comparison  purposes,  the  transfer  functions  for  the  three  systems  are 
given  in  Table  2.  The  numerator  coefficients  and  are  functions 

of  vehicle  and  sight  characteristics,  and  therefore  depend  on  tracking  range 
or,  equivalently,  Ty,  On  the  other  hand,  except  for  the  sight  time  constant 
Ty  in  the  typical-sight  case,  the  denominator  is  a  function  of  vehicle  charac¬ 
teristics  alone. 

The  optimal -control  model  (0CM)  of  the  pilo£5alssumes  that  the  well-trained, 
well -motivated  human  operator  chooses  his  control  inputs  up,  subject  to  human 
limitations,  such  that  following  objective  function  is  minimized 

cT 


1  -  c  r  1  ^  1 

Jp  '  E{T-~  T 


(y'Qy  +  UpRUp  +  UpGup)dt} 


(3) 


and  G  is  selected  to  obtain  a  chosen  neuromuscular  lag  time  constant  r^. 
pilot’s  input  is  then  expressed  in  the  scalar  case  as 

t  kjU  =  -  K  x  -  un 
N  p  x  p 


The 


L 


n 


The  selected  parameters  of  the  model  in  this  investigation  are  given  in  Table 
3  ,  and  are  consistent  with  those  in  Harvey's  investigation  of  this  task. 
(Note  that  nominal  values  of  model  parameters  have  been  chosen  which  do  not 
reflect  the  pilot's  acceleration  environment.  This  was  done  to  be  consistent 
with  fixed-base  simulation  results.) 

TABLE  2 

_ TRANSFER  FUNCTION  COMPARISON _ 

Second  order  plant: 

dJL  ,  1LZ 

s2 

Ideal  sight  display: 

cjs)  ,  *(",s3  »  n2s2  ♦  n3s  *  v _ 

Typical  si  gilt  display: 

cjs)  ,  K(H,s3  *  m/  ♦  M35  ♦  M4) _  __ 

s2(Tfs  *1)<V  *  l)(s2  7  2SAps  77[p) 

K  =  -16  t  /0 
6 

",  -  -962Vf 

Vif'WW 

n3  ■  V", *  +  <°  *  V'V.  -  w 

N4  ■  "„ZS  -  Za"« 

",  *  -  -962! 

"2  ■  N2  ‘  Tf -  ■962><V"«  -  V*> 


M3  =  N3;  M4  =  N4 


TABLE  3 

PILOT  MODEL  PARAMETERS 


Observation  delay,  t  =  0.2  sec 
Neuromuscular  time  constant,  =  0.1  sec 

Observed  variables,  y'  =  [e.e.a.a,] 

Cost  function  weightings,  Q  =  [16,1,0,4] 

yii 

R  =  0 
u 

Full  attentional  allocation 

Observation  thresholds,  T  =  T  =  0.65  deg 

F  A  3 

T*  =  T^  =  1.3  deg/sec 
2 

Observation  noise  variance  V  =  up  o  ,  p  =  .01 

yi  y  y,  y 
2 

Motor  noise  variance,  V  =  up  o  .  ,  p  =  .001 

up  u  up  u 

Shown  in  Table  4  is  the  predicted  (ms)  performance  using  this  0CM  as 
compared  with  Harvey's  (fixed-base  simulator)  experimental  results  for  the 
typical  sight  display.  Also  shown  is  the  predicted  performance  for  the 
idealized  sight  for  reference. 

SYSTEM  AUGMENTATION 

With  the  above  system  models,  we  represent  the  system,  including  the 
pilot's  control  input  as  well  as  the  augmentation  input,  by  the  relation 

x  =  Ax  +  Bu  +  Bu  +  w 
p  aug 

where  u  is  an  equivalent  stick  input,  a  scalar,  in  this  case.  As  shown 
o  uy 

previouslyP^if  u_  is  chosen  to  minimize  the  objective  function 

a  Uy 


Mil 


=  1000  ft  a 


Experiment 

Analytical 
(Typical  sight) 

Analytical 
(Ideal  sight) 


=  1000  ft  0j  =  5  g 


Experiment 

Analytical 
(Typical  sight) 

Analytical 
(Ideal  sight) 


3000  ft  o 


Experiment 

Analytical 
(Typical  sight) 

Analytical 
(Ideal  sight) 


3000  ft  a 


Experiment 

Analytical 
(Typical  sight) 

Analytical 
(Ideal  sight) 


TABLE  4 

PERFORMANCE  COMPARISON 


Parameter  (rms) 
£(deg)  A(deg) 


q(deg/sec) 


6r(deg) 


i  8.0 

2.1 

i  8.7 

2.1 

i  7.9 

2.0 

2.4 

8.7 

6.4 

1.8 

2.3 

10.1 

7.3 

1.8 

1.1 

9.7 

6.0 

1  .6 

aug  p 


J  +  EU""  t  I  u'  Fu  dt> 

T-m=  [  J 


aug  aug 


where  Jp  is  given  in  Eqn.  3  ,  the  relation  for  the  augmentation  input  is 

u  =  -F~Vk.  x  -  F'Vk.  u 

rt2 


(The  matrices  and  are  obtained,  as  discussed  in  Reference  1,  by 
simultaneously  solving  two  coupled  Riccati  equations,  one  yielding  the 
augmentation  gains,  and  one  yielding  the  pilot's  gains.). 

The  result,  of  course,  is  the  pilot  controlling  the  augmented  system, 
which  is  now  described  by 

x  =  (A  -  BKjx  +  B( I  -  Ku)up  +  w 
or 

x  =  A  x  +  B  u  +  w 
P  P  P 


Second-Order  Plant 

Let  us  initially  consider  the  effects  of  pilot  optimal  augmentation  of 
the  second  order  plant  (for  which  we  may  define  x'  =  (@c,  §c,  e,  e)  and 
y'  =  (e,  e)).  Shown  in  Fig.  2  is  the  augmented  open- loop  and  pilot  closed- 
loop  root  locus  obtained  via  optimal  augmentation  for  various  uau^  objective- 
function  weights,  or  F  in  Eqn.  4  .  (A-E  indicate  open-  and  closed-loop  cases 

with  increasing  augmentation.)  Note  the  shape  of  the  augmented  (open-loop) 
system  root  locus  as  well  as  that  of  the  pilot  closed-loop  system  corresponds 
to  the  know£6c)symptotic  behavior  of  the  closed-loop  roots  of  a  linear  optimally 
controlled  single-input,  single-output  system.  In  this  case,  the  root  locus 
assumes  a  Butterworth  configuration  of  order  two.  (This  results  from  two  open- 
loop  plant  poles  and  no  zeros.  Except  the  pilot  closed-loop  poles  include  a 
third  introduced  via  the  pilots  neuromuscular  lag.) 


ocus-Second  Order  Plant 


w 


The  distance  between  open-loop  and  closed-loop  eigenvalues  is  a  maxi¬ 
mum  for  no  augmentation  (i.e.  while  the  open-loop  and  closed-loop  poles 

coalesce  for  increasing  augmentation,  eventually  eliminating  the  pilot  compen¬ 
sation,  and  the  system  monotoically  approaches  a  pure  automatically  controlled 
system. 

To  aid  in  further  interpretation  of  these  results,  consider  the  human 

operator  describing  function  in  a  compensatory  task  to  be  approximated  by 
K  e‘xS(T,s  +  1) 

H(s)  ■  p(V'o 

or  a  gain,  low-frequency  lead,  effective  delay,  and  neuromuscular  lag.  Our 
augmented  or  unaugmented  second-order  plants  may  be  described  in  general  by  the 
transfer  function 

K 


G(s)  = 


__ 

s  +  2^'jS  +  u> 


The  piloted  closed-loop  denominator  now  becomes  (ignorinq  the  delay) 

s3  +  (—  +  2cw)s2  +  [oj2  +  —  (2e&»  +  K K  T.  )]s  +  —  {J  +  KK  ) 
tN  tN  p  L  tn  p 

If  we  write  this  closed-loop  denominator  as 
(Tcs  +  1 )  ( s2  +  2?cucs  +  co2 ) 


or 

S3  +  (i-  +  2c  u>  )s2  +  (“r  +  w  )s  +  w2 

'c  c  u  c  c  c  c  c  C 

we  see  that  the  pilot's  gain  Kp  is  related  to  the  difference  between  open-  and 

closed-loop  frequencies  by  the  relation 


_  TN  2 

-  r  “c 


c 


2 

(0 
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Therefore,  the  pilot's  gain  tends  to  be  proportional  to  the  difference  be¬ 
tween  the  squares  of  closed-  and  open-loop  natural  frequencies ,  or  equiva¬ 
lently,  the  squares  of  distances  of  the  poles  from  the  origin. 


Also,  the  pilot's  lead  T^,  is  related  to  the  open-  and  closed-loop 
parameters  by  the  relation 

T  =  tN[(“c  -  +  2<WTc  -  ^ 

L  ,tN  2  2 , 

It-  “  ) 

C 


Therefore,  the  pilot's  lead  would  tend  to  increase  with  increasing  a(cw) 

2 

/Aw  between  open-  and  closed-loop  eigenvalues.  With  these  observations,  we 
are  able  to  qual i tiati vely  relate  the  trends  from  the  root  locus  (Fig.  2)  to 
the  parameters  in  the  classical  human  operator  describing  function. 

In  the  case  of  higher  order  systems,  however,  the  pilot  describing 
function  may  not  be  well  approximated  by  the  simpler  forms,  and  the  relation¬ 
ships  between  the  root  locus  characteristics  and  the  predicted  human  describing 
function  parameters  are  much  more  ^omplex  and  difficult  to  express  quantita¬ 
tively.  Still,  the  effects  of  augmentation,  and  of  the  pilot's  gains  are  de¬ 
picted  in  the  root  loci,  and  the  difference  between  the  open-  and  closed-loop 
roots  indicate  at  least  qualitatively  the  required  level  of  human  equalization 
involved.  And  finally,  the  augmented  (open-loop)  system  root  locus  yields  im¬ 
portant  information  about  the  desirable  trends  in  the  systems  dynamics. 


eigenvalues  are  shown  in  Fig.  3  . 


The  open-loop  roots  include,  the  unaugmented  case,  the  two  roots  at  the 
origin,  the  vehicle  short-period  roots,  and  the  actuator  root  at  -20(rad/sec. ) 

For  the  piloted  closed-loop  system,  still  unaugmented,  the  roots  include 
a  low-frequency  pair,  a  high-frequency  pair,  and  two  real  roots,  one  at 
-3.7  rad/sec  and  one  (not  shown  for  clarity)  at  -20  rad/sec. 

As  the  level  of  augmentation  is  increased,  as  indicated  by  cases  A-C, 
the  frequency  and  damping  of  the  vehicle  short  period  mode  is  increased,  and 
the  two  poles  originally  at  the  origin  move  toward  the  piloted  closed-loop 
low-frequency  pair  due  to  the  closure  of  the  tracking  loop  (i.e.,  feedback 
of  e  and  c).  With  the  pilot  closing  the  loop  on  the  augmented  systems  a 
unique  pair  of  high  frequency  roots  are  associated  with  each  level  of  augmen¬ 
tation,  and  this  pair  appears  to  follow  a  Butterworth  pattern.  However,  for 
the  levels  of  augmentation  considered,  the  two  real  closed-loop  roots  (at 
approximately  -20  and  -3.7  rad/sec)  are  relatively  insensitive  to  augmenta¬ 
tion  level.  Also,  the  low-frequency  pair  (at  about  -.7  +  .7j)  are  invariant 
with  augmentation.  This  pair  also  corresponds  closely  to  a  pair  of  plant 
zeros . 

This  behavior  is  consistent  with  the  asymptotic  properties  of  linear 
optimally  controlled  system  roots  cited  previously.  For  single-input,  single¬ 
output  systems  p  of  the  roots  tend  toward  the  p  open-loop  zeros  (in  the  left 
half  plane),  and  the  remaining  roots  assume  a  Butterworth  configuration  of 
order  n-p.  Furthermore,  the  systems  with  input  and  output  vectors  of  equal 
dimension  (i.e.,  systems  with  square  transfer-function  matrices),  p  of  the 
closed-loop  roots  tend  to  the  p  "zeros"  of  the  determinant  of  the  transfer- 
function  matrix.  Although,  the  transfer-function  matrix  is  not  square  in  this 
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case,  the  invariance  of  any  closed-loop  roots  may  still  be  explained  in 
terms  of  the  above  properties  and  the  existence  of  system  “pseudo  zeros." 

Shown  in  Fig.  4  is  the  comparison  of  these  (vehicle  short  period) 
results  with  those  of  Hollis.  In  his  analysis,  Hollis  determined  the  ve¬ 
hicle  augmentation  yielding  the  best  pilot  rating  via  the  "paper-pilot". 

Pitch  rate  and  plunge  acceleration  were  used  for  augmenting  only  the 
vehicle  with  a  fixed  sight  (with  no  display  dynamics).  It  is  seen  that  the 
results  agree  with  the  trend  determined  by  Hollis'  consideration  of  several 
sets  of  vehicle  dynamics.  In  all  cases,  the  "optimal"  augmentation  in¬ 
creased  the  frequency  and  damping  of  this  vehicle  mode.  This  trend,  as 
noted  by  Hollis,  also  agrees  with  the  military  specifications  on  flying 
qualities.  (As  we  shall  see  later,  this  is  not  the  case  for  the  other  sight 
display. ) 

The  predicted  improvement  in  (rms)  trackinq  accuracy  for  the  ideal 
sight  is  shown  in  Table  5,  along  with  the  rms  statistics  on  stick  deflection 
and  rate  of  deflection  (i.e.,  physical  workload).  Recall  that  the  stick 
deflection  was  not  considered  penalized  in  the  pilot's  objective  function, 
hence  little  reduction  in  this  parameter  is  obtained  through  augmentation. 
However,  rms  stick  rate  is  significantly  reduced. 

Finally,  Table  6  lists  the  optimum  full-state  feedback  gains  in  this 
idealized  case,  for  each  level  of  augmentation. 

For  this  sight  display,  increasing  the  tracking  range  D,  from  1000  ft. 
to  3000  ft.  changed  the  results  very  slightly.  The  feedback  gains  varied 
approximately  15  percent  and  the  augmented  and  piloted  system  root  loci  were 
almost  identical  to  that  in  Fig.  3.  The  position  of  the  invariant,  low-frequency 
piloted-system  roots  did  change  from  -.7  +  . 7j  to  -.3  and  -1.2  rad/sec,  but 
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Vehicle 


the  augmentation  moved  the  poles  at  the  origin  toward  those  roots  as  in  the 
1000  ft.  case. 


Typical  Sight  Display 

Considering  now  the  system  with  the  more  typical  display  dynamics  (lead 
angle  computed  with  Eqn.  2  ),  the  open-  and  closed-loop  root  loci  are 
shown  in  Fig.  5.  The  unaugmented  system  roots  in  this  case  include  those 
considered  above,  plus  the  lead-angle  time  constant  T,.  Comparing  these  re¬ 
sults  with  those  in  Fig.  3,  we  see  an  entirely  different  trend  in  the  vehicle 
short-period  roots!  With  the  more  sensitivity,  or  "active,"  sight  display 
in  this  case,  the  optimum  short  period  frequency  remains  relatively  con¬ 
stant  while  the  damping  is  increased.  Recall  that  the  results  obtained 
by  Hollis  included  no  display  dynamics,  and  that  the  results  agreed  with 
those  obtained  with  the  less  active,  ideal  sight  display.  Additionally, 
the  augmented  short-period  roots  in  the  case  considered  here  (Fig.  5) 
approach  an  invariant  pair  of  piloted  closed-loop  eigenvalues,  a  trend  not 
observed  in  the  ideal-sight  case. 

The  behavior  of  the  roots  near  the  origin  is  as  observed  previously, 
the  augmentation  closing  the  tracking  loop  and  moving  the  plant  poles  at 
the  origin  toward  the  two  invariant  closed-loop  eigenvalues.  Also,  the 
actuator  and  sight  time  constants  move  together,  due  to  augmentation,  and 
eventually  become  a  high-frequency,  complex  pair.  Finally,  as  before,  a 
unique,  high-frequency  closed-loop  pair  is  associated  with  each  level  of 
augmentation,  and  approaches  a  Butterworth  pattern.  A  seventh  closed-loop 
root  at  -20  rad/sec  is  the  actuator  root  remaining,  but  is  not  shown  for 
clarity. 
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Root  Locus-Typical  Sight  (1000  ft) 


Increasing  the  range  to  the  target  D  from  1000,  to  3000  ft.  results  in 
the  root  loci  shown  in  Fig.  6.  Now,  not  only  do  the  plant  zeros  (and  there¬ 
fore  the  "pseudo-zeros")  depend  on  D,  but  the  sight  time  constant  is 
affected  as  well.  As  a  result,  this  root  locus  is  appreciably  different 
from  that  in  Fig.  5,  particularly  at  the  higher  level  of  augmentation. 

Still  present,  however,  are  the  four  invariant  closed-loop  poles,  two  real 
and  two  complex  in  this  case,  and  the  family  of  complex  closed-loop  roots 
at  the  higher  frequency.  Note  in  this  case,  the  migration  of  the  sight  time 
constant  with  augmentation  toward  the  closed-loop  at  -.3  rad/sec,  and  the 
behavior  of  the  augmented  short-period  roots  for  level  C. 

The  effect  of  augmentation  (level  B,  D  =  1000  ft)  on  the  system's 
frequency  response  (e/6)  is  depicted  in  Fig.  7.  Closure  of  the  tracking 
loop  (via  e  and  e  feedback)  eliminates  the  1/s  characteristic,  at  low  fre¬ 
quency,  while  the  increased  damping  of  the  short  period  mode  (at  w  *  3.  rad / 
sec)  is  evident.  The  effect  of  this  level  of  augmentation  on  the  predicted 
frequency  response  of  the  pilot  is  shown  in  Fig.  8.  With  augmentation,  the 
pilot  appears  to  increase  his  gain,  as  well  as  adding  slightly  more  low 
frequency  lag.  Also,  his  required  lead  near  u  =  3  rad/sec.  would  appear  to 
be  decreased. 

The  pilot  frequency  response  shown  in  this  figure  is  that  of  an 
equivalent  describing  function  formed  from  the  transfer  function  matrix 
(H(s))  obtained  from  the  optimal  control  pilot  model.  This  equivalent  de¬ 
scribing  function  is  analogous  to  the  classical  human  describing  function 
measured  experimentally  in  a  compensatory  task,  and  is  obtained  from  the 
relation 


«st(s) 


rrsr 


Equiv. 


sst(s) 

7TTT 


+  s 


6st(s) 

c(5) 


=  hn(s) 


sh12(s) 


26 


Root  Locus-Typical  Sight  (3000  ft) 


Augmented 
(level  B) 


and  where  6$t(s)  =  H(s)  y(s). 

The  effect  of  the  augmentation  on  the  predicted  rms  performance  is 
shown  in  Table  7.  The  rms  error  e,  stick  deflection  and  stick  rate  are 
given  for  D  =  1000  ft.,  oaj  =  5  g's.  Here  again  we  see  a  significant  im¬ 
provement  in  tracking  error  while  the  rms  stick  activity  remains  approximately 
constant  with  increased  augmentation.  As  before,  however,  we  note  that 
stick  activity  was  not  penalized  in  the  objective  functions.  The  pilot  stick 
rate  is  also  reduced  with  augmentation. 

The  above  improvements  in  tracking  accuracy  were  obtained  while  the 
predicted  elevator  activity  actually  decreased  slightly  with  increased 
augmentation  (2.96  deg-rms  unaugmented  to  2.81  deg-rms  for  level  C).  This 
is  due  to  the  fact  that  the  augmentation  includes  pilot  stick  position 

feedback  (K.  ),  so  the  net  commanded  elevator  deflection  by  the  pilot  be- 

6st 

comes 


«  K(1  -  K  ) 
E  6st 


St 


where  K  is  the  actual  stick  gain.  With  his  effective  stick  gain  reduced, 
the  pilot's  contribution  to  commanded  elevator  is  reduced  while  the 
augmentation's  contribution  is  increased.  The  sum  of  the  two  results  in  the 
relatively  constant  net  rms  elevator  with  increasing  augmentation.  The 
optimum  gains  for  the  vehicle  states  for  the  three  levels  of  augmentation 
is  given  in  Table  8. 
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SECTION  HI 


MULTI-AXIS  ANALYSIS 
AND  AUGMENTATION 

The  more  challenging  (and  more  realistic)  problem  of  pilot-vehicle 
analysis  and  augmentation  for  a  multi -axis  air-to-air  tracking  task  will 
now  be  considered.  The  dynamics  of  the  system,  as  in  the  pitch  tracking 
analysis  in  the  previous  section,  include  both  vehicle  and  display 
dynamics.  However,  the  flight  condition  to  be  considered  is  a  highly 
banked  turn  with  a  normal  load  factor  of  four  (i.e.,  4  g  turn).  This 
involves  several  issues  not  frequently,  and  in  some  cases  never,  considered 
in  previous  investigations  of  pilot-vehicle  dynamics.  These  extenuating 
issues  result  from  the  significant  amount  of  unsymmetrical  coupling  between 
the  elevation  and  azimuth  axes  of  the  system  and  the  multiple  control  inputs 
involved  (i.e.,  elevator,  aileron,  and  rudder). 

Previous  investigations  of  multi-axis  tasks  have  been  few  in  number 
and  those  that  have  considered  the  effects  of  multi-axis  aircraft  control^8”10-^ 
have  actually  treated  the  dynamics  of  both  axes  independently,  reducing  the 
man's  attentional  allocation  to  each  axis  due  to  the  simultaneous  control 
task.  When  applicable,  this  approach  is  nice  since  the  dynamic  model  is 
smaller  and  more  manageable.  And  in  wings-level  flight  in  a  position 
tracking  task  the  dynamics  are  uncoupled,  hence  independent.  However,  in 
a  highly-banked  turn,  and  especially  in  an  anqular tracking  or  pointing  task, 
the  axes  are  highly  interactive. 

The  problems  to  be  addressed  specifically  in  this  section  then  include 
1)  the  determination  of  any  necessary  chang°s  in  the  pilot  model  parameters 
for  multi-axis  modeling  from  those  determined  from  the  pitch  tracking 
analysis  of  Harvey^lnd  2)  the  application  of  the  pilot-optimal  control 
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synthesis  procedure  to  this  multi-axis  analysis  and  prediction  of  the 
improvement  in  performance,  if  any. 

To  accomplish  the  first  objective,  the  comparison  will  be  made  betv'en 
analytical,  or  model,  results  and  simulation  data.  The  data  was  obtained 
previously  in  the  large-amplitude,  motion  simulator  of  the  Flight  Dynamics 
Laboratory,  Wright-Patterson  Air  Force  Base  as  part  of  another  study  performed 
by  the  Flight  Dynamics  Laboratory.  The  simulation  involved  F106  aircraft 
dynamics,  and  investigated  several  sets  of  display,  or  sight  dynamics.  The 
display  to  be  considered  here  is  the  ideal  display  (or  perfect  director 
sight)  considered  in  the  previous  section.  The  development  of  linearized 
mathematical  models  for  the  vehicle  and  display  in  this  multi-axis  configura¬ 
tion  is  presented  in  Appendices  B  and  D  of  this  report. 

The  approach  to  the  multi-axis  pilot  modeling  task  was  to  attempt  to 
match  analytical  and  simulated  statistical  results  in  the  form  of  rms  values 
of  system  states,  outputs,  and  control  inputs.  The  initial  objective  function 
was  a  direct  extension  of  the  cost  function  determined  experimentally  by 
Harvey  -  that  used  in  the  previous  section.  The  observation  errors,  neuromotor 
dynamics  and  number  of  control  inputs  were  then  varied  to  obtain  the  best 
match  (via  eyeball  fit)  of  rms  values.  As  will  be  seen,  some  additional 
parameters  needed  variation  and  describinn  functions  were  reauired  in  some  cases. 

Significant  to  this  discussion  of  pilot  modeling  is  the  fact  that  the 
simulated  task  involved  tracking  a  non-maneuvering  tarqet.  The  pilot  rolled 
into  the  turn, attempted  to  capture  the  target  in  his  display,  and  then 
tracked  in  approximately  a  steady-state  condition.  In  terms  of  a  linear 
model,  this  corresponds  to  a  system  with  no  process  driving  noise.  In  fact, 
the  only  contribution  to  tracking  error  in  this  case  is  the  pilots  observation 
and  neuromotor  noise.  To  attempt  to  model  this  situation,  additive  motor 


noise  of  specified  variance  was  injected  into  the  pilot's  control  input, 
and  the  variance  of  this  noise  was  chosen  so  as  to  match  model  and  simulated 
rms  control  activity.  The  observation  noise  statistics  were  treated  in 
the  usual  manner,  that  is,  the  variance  of  observation  noise  was  specified 
in  terms  of  a  noise- to-signal  ratio  (e.g.  -20  dB  for  full  attentional 
allocation  to  a  particular  output  quantity). 

SYSTEM  MODEL 

For  this  "frozen-point"  analysis,  the  system  dynamics  are  linearized 
about  a  four-g,  turning  flight  condition,  the  mathematical  definition  of 
which  is  given  below. 

TABLE  9 

STEADY  STATE  PARAMETERS 

Altitude,  h  -  10,000  ft;  Mach  =  0.72 

Flight  Velocity,  Uss  =  775  ft/sec 

Flight  path  angle,  Oss  =  0. 

Bank  angle,  =  75.5  deg 

Normal  acceleration,  Az  =  -4g  =  -  128.8  ft/sec 

Turn  rate,  v  =  0.161  rad/sec 

Pitch  rate,  Q$s  =  0.156  rad/sec 

Yaw  rate,  R$s  =  0.040  rad/sec 

Angle  of  attack,  a$s  =  14.5  deg 

Elevator  deflection,  a  =  -20.2  deg 

tss 


Furthermore,  unlike  a  level  flight  condition,  the  engagement  geometry 
requires  a  non-zero  steady-state  line  of  sight  to  the  target,  sight  dis¬ 
played  lead  angle,  and  relative  heading  between  attacker  and  target.  Since 
the  steady-state  line  of  sight  must  be  a  constant  for  this  (statistically 
stationary)  modeling  approach,  the  engagement  geometry  is  essentially 
specified  once  the  tracking  range  is  selected.  The  resulting  steady- 
state  values  may  then  be  compared  to  the  average  values  obtained  from 
simulation  to  verify  that  the  simulation  actually  reflected  a  steady-state 
tracking  condition.  Finally,  note  that  the  steady-state  tracking  error 
must  be  zero  for  the  perturbation  quantities  obtained  from  the  linear 
mathematical  model  to  be  representative  of  total  tracking  error. 

From  the  1 ines-of-sight  equations  developed  in  Appendix  A  and  the 
display  lead  angle  equation  in  Appendix  B,  the  engagement  geometry  may 
be  specified  as  in  Table  10  for  both  aircraft  on  a  common  circular  trajectory. 

The  steady-state  lead  angle  may  be  found  from  the  relations  presented 
in  Appendix  B.  Assuming  a  projectile  muzzle  velocity  of  3300  ft/sec,  a 
time  of  flight  of  0.786  sec.,  and  a  ballistic  jump  correction,  Jv,  of  0.185, 
yields  the  following  lead  angles  in  the  steady-state  flight  condition. 


'El 


ss 


AZ 


0.13  rad 


-0.03  rad 


ss 


These  lead  angles  are  as  shown  in  Figure  9  alonn  with  the  definition  of 
the  components  of  the  line  of  sight  to  the  target.  Finally,  the  steady-state 
tracking  error  is  taken  to  be 

£nss  =  “elss  -  xeiss  ‘ 0 


=  0 
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TABLE  10 


ENGAGEMENT  PARAMETERS 


Target  altitude,  h-j.  =  10,000  ft;  Mach  =  0.72 

Target  Velocity,  U  =  775  ft/sec 
s  T 

Target  flight  path  angle,  0-^  =  0 
Target  bank  angle,  =  75.5  deg 
Distance  to  target,  d  =  2000  ft 
Target/Attacker  relative  heading,  ay  =  23.6° 

2 

Target  normal  acceleration,  A7  =  -128.8  ft/sec 

lT 

Target  acceleration  in  attackers  coordinates: 

At  =  -50.  ft/sec2 

'x 

At  =  -3.  ft/s ec2 
'y 

At  =  -119.  ft/sec2 

'z 

Target  velocity  vector  in  attacker's  coordinates: 

VT  =710  ft/sec 

'x 

VT  =  78  ft/sec 
'Y 

VT  =  -300  ft/sec 

'z 


The  comparison  of  the  above  "theoretical”  steady-state  with  the  mean 
value  of  the  parameters  obtained  from  the  two  simulation  runs  are  presented 
in  Table  11.  (These  two  runs  were  the  only  cases  for  the  modeled  flight 
condition,  display  dynamics,  etc.)  It  is  apparent  that  Run  109  did  not 
simulate  the  modeled  steady-state  situation,  and  checking  the  actual  time 
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TABLE  11 

STEADY-STATE  PARAMETER 
COMPARISON 


- 

Tracking  Errors (rad) 

Lead  Angles  (rad) 

£ELss 

£AZ  c 
ss 

*E’ss 

XAZ 

ss 

Theoretical 

0.0 

.127 

-.031 

Run  133 

.026 

.123 

.021 

.025 

-.028 

Target  Velocity 

Components  (gun  coord.) 

2 

Attacker  Acceleration  (ft/sec  ) 
(body  axis) 

\s 

Wy  ft/sec 
ss 

Ax 

Ay 

Az 

Theoretical 

78. 

-300. 

32. 

0.  - 

125. 

Run  133 

101. 

-312. 

12. 

-1.  - 

104. 

Run  109 

27. 

-  65. 

8. 

-3.  - 

42. 

Body  Euler  Angles  (deg.) 

Body  Angular  Rates  (deg/sec) 

°ss 

<t> 

ss 

P 

ss 

“ss 

R 

ss 

Theoretical 

3.6 

75.5 

-  0.6 

8.9 

2.2 

Run  133 

3.6 

72.6 

-  0.6 

8.3 

2.5 

Run  109 

3.4 

28.9 

38.0 

2.9 

2.6 

Velocity  Angles  (deg. ) 

Control 

Deflections  (in.) 

ass 

i 

6ss 

6e 

st 

6a  f 

st 

^0 

Kped 

Theoretical 

14.5 

0.0 

-1. 

0. 

0. 

Run  133 

10.9 

0.6 

-2.3 

-0.4 

0.2 

Run  109 

3.8 

0.6 

-1.0 

-0.1 

1.3 
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histories  of  the  parameters  verified  that  a  steady-state  situation  was 
in  fact  not  simulated.  That  is,  rather  than  constant  mean  values,  the 
means  were  time  varying  for  all  parameters  and  a  statistically  stationary 
model  is  therefore  invalid. 

On  the  other  hand,  the  data  from  Run  133  appears  to  agree  quite  well. 
The  only  discrepancy  is  the  difference  in  range  to  the  target  (not  given 
in  the  table).  The  assumed  encounter  geometry  considered  a  constant  range 
of  2000  feet,  while  the  simulated  range  varied  monotonically  from  1800  to 
3000  ft,  with  a  mean  value  of  2315  feet.  The  primary  effect  of  this 
variation  is  on  the  standard  deviation  of  the  displayed  lead  angle  about 
the  mean,  the  standard  deviation  is  much  larger  than  would  be  the  case 
with  constant  range.  With  this  in  mind,  however.  Run  133  will  be  taken  as 
the  simulation  results  to  be  compared  with  the  analytical  statistics. 

The  linear  dynamic  model  may  now  be  defined,  and  expressed  in  the 
familiar  form 

x  =  Ax  +  Bup 

with  the  state  vector  of  perturbation  variables  taken  as 

x  —  Vj,  Wy ,  d,  ^9  ^ 9  & 9  ^9  ^ J 

where 

=  relative  heading  angle  Uy  - 
Vy,  Wy  =  target  velocity  in  attacker's  coordinates 
d  =  range  variation 

eEl’  eAZ  =  ^ne  angle  components 

e,  $  =  attacker  Euler  angles 


a,  6  =  angles  of  attack  and  sideslip 


p,  q,  r  =  stability  axis  angular  rates. 


The  pilot's  control  vector  is 


Up  =  «  <$»  ,  -I 

st  Ast  RpedJ 


where  as  with  the  states,  the  stick  and  pedal  deflections  are  perturbations 
from  their  steady  state  values. 

The  pilot's  observation  vector  was  ultimately  taken  to  be 


with 


where 


yp ( t  +  t)  =  Cx(t)  +  Dup(t) 


— T  r  •  •  •  •  -  _ 

yp  »  UEv  eEl>  eAZ’  eAZ*  XE1’  XEl  *  XAZ’  XAZ’  Ay’ 


cEl  ’  eAZ  =  e1evat‘ion  and  azimuth  tracking  error 
XEl  ’  XAZ  =  e^eva^on  and  azimuth  displayed  lead  angle 
Ay  =  lateral  acceleration  at  the  cockpit 
4>  =  bank  angle  (or  relative  bank  angle) 


Note  that  the  errors  and  lead  angles  are  included  as  in  the  pitch  tracking 
analysis  of  the  previous  section.  Lateral  acceleration  and  bank  angle 
must  now  also  be  included  in  the  multi-axis  task  with  motion.  These  latter 
variables  didn't  need  consideration  in  the  previous  analysis  as  it  was 
based  on  a  fixed-base  pitch  tracking  experiment.  They  definitely  must 
be  included,  however,  in  this  more  complex  task  with  motion  simulation  as 
we  shall  see  later. 

With  the  steady-state  parameters  developed  previously,  and  the 
perturbation  equations  developed  in  Appendices  A,  B,  C,  and  0,  the  matrices 
A,  B,  C,  and  0  in  the  above  equation  are  as  given  in  Appendix  E. 


To  be  addressed  now  is  the  coupling  between  the  longitudinal  and 
lateral  directional  axes,  and  the  piloting  "strategy"  in  this  multi-axis 
task.  Shown  in  Figure  9  previously  was  the  schematic  of  the  pilot's 
view  through  the  heads-up  display.  Now  one  classical  approach  is  to  null 
the  azimuth  error  by  rolling  the  aircraft  (using  aileron  and  rudder)  such 
that  all  the  error  is  in  the  pitch  axis. 

This  situation  is  depicted  in  Figure  in  ,  where  for  simplicity  of 
explanation,  the  sight  (pipper  and  retical)  is  assumed  fixed  in  space  and 
the  pilot  has  rolled  about  the  line  in  space  represented  by  the  sight 
position.  It  is  apparent  here  that  by  using  lateral-directional  control, 
the  elevation,  or  pitch  error  has  been  increased.  Conversely,  a  purely 
pitching  maneuver  introduces  no  change  in  azimuth  error  (assuming  the  sight 
to  be  fixed). 

This  anti -symmetric  cross-coupling  is  depicted  graphically  in  Figure 
11  ,  which  shows  the  frequency  response  of  the  system's  anti -symmetric 
transfer  functions.  The  coupling  between  aileron  input  and  elevation 
error  is  20  to  40  dB  larger  than  that  between  elevator  input  and  azimuth 
error.  Here,  of  course,  the  complete  dynamics  of  the  vehicle  and  display 
are  included,  while  in  the  previous  discussion  of  the  display  schematics 
we  assumed  the  sight  fixed  in  the  display.  System  eigenvalues  and  frequency 
responses  of  the  remaining  transfer  functions  will  be  discussed  later 
with  augmentation  results. 

PILOT-VEHICLE  ANALYSIS 

The  baseline  pilot  model  was  a  direct  extension  of  the  pitch  tracking 
analysis  presented  previously,  and  this  model  is  summarized  in  Table  12. 
Reiterating  the  approach  outlined  previously,  the  cost  function  weightings 
in  both  axes  are  idential  to  those  for  the  pitch-tracking  case.  The 
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ect  of  Rol 1 ing 


on  Geometry 


BASELINE  PILOT  MODEL 


Observation  Vector,  y  C^ei  *  CE1*  ^Az*  eAz*  ^El*  ^Az*  ^Az^ 

Objective  Function  Weights,  Q  =  16. ,  Q*  =  1 . ,  Q*  =  4. 

(both  axes) 

Observation  Thresholds,  T  =  T  =  0.05  deg,  T*  =  T?  =  0.10  deg/sec 

Observation  Noise  Ratio,  =-20  dB  for  full  attention 

Fractional  Attention,  f.  =  0.5  all  observed  variables 

Observation  delay,  t  =  0.2  sec 

Neuromuscular  lag,  =  0.2  sec  all  inputs 

Motor  Noise  Variance,  varied  to  match  rms  control  inputs 

Control  inputs,  ,  6.  only 

stick  stick 

motor  noise  variance  was  modified  to  match  model  and  simulated  rms  control 
activity,  and  the  parameters  determining  observation  noise  (i.e.  observation 
thresholds  and  attentional  allocation)  along  with  neuromuscular  lag  were 
varied  in  the  attempt  to  match  the  statistics  of  the  system  states  and 
outputs. 

Initially,  only  two  control  inputs  were  considered,  elevator  and 
aileron  stick  deflection,  as  the  importance  of  the  rudder  in  this  task 
was  not  clear  a  priori.  Part  of  the  modeling  effort  then  was  to  determine 
the  necessity  for  including  the  rudder,  and  it  will  be  shown  that  inclusion 
is  warranted. 

Finally,  recall  that  the  target  is  not  maneuvering,  so  the  magnitude 


of  the  motor  noise  variances  must  be  specified  rather  than  noise-to- 
signal  ratio,  the  more  typical  case.  The  resulting  noise-to-signal 
ratio  obtained  from  matching  control  activity  is  then  one  of  the  modeling 
results  in  this  analysis. 

The  comparison  of  the  baseline  model  results  with  those  obtained 
from  simulation  is  shown  in  Table  13  .  Considering  the  use  of  essentially 
a  "single  axis"  model,  the  results  are  amazingly  good,  the  largest 
differences  appearing  in  lead  angles,  bank  (<j>)  and  side-slip  angle  (s), 
and  higher  angular  rates  resulting  from  the  model. 

To  assess  the  first-order  sensitivities  of  the  rms  performance  to 

increased  observation  errors  and  pilot  bandwidth  (i.e.  neuromuscular 

lag),  two  "perturbations"  to  the  baseline  model  were  evaluated.  The 

first  involved  reducing  the  fractional  attention  allocation  from  0.5,  the 

inverse  of  the  number  of  controlled  axes,  to  0.125,  the  inverse  of  the 

dimension  of  the  observation  vector.  Additionally,  the  threshold  on 

angular  rate  was  increased  from  0.1  to  0.18  deg/sec  consistent  with 

[11] 

some  earlier  studies.  The  second  case  considered  the  above  changes  plus 
increased  the  neuromuscular  lag  time  constant  from  0.2  seconds  to  0.3 
seconds  for  both  control  inputs. 

The  results  for  these  two  additional  models  are  shown  in  Table  14, 
along  with  the  baseline  model  results.  The  effect  of  the  increase  in  the 
observation  noise  due  to  reduced  attentional  capacity  and  increased 
threshold  is  an  increase  in  all  rms  values  of  states  and  outputs.  But 
this  is  as  expected  since  increased  observation  noise  is  similar  to 
increased  motor  noise  in  effect.  Note,  however,  that  the  increases  appear 
to  be  greater  in  the  lateral  axis,  with  rms  azimuth  error  (e^  )  increasing 
65  percent  and  side-slip  angle  (b)  and  roll  rate  (p)  increase  10  percent. 


45 


It  would  appear  that  the  lateral-directional  axis,  then,  is  more  sensitive 
to  observation  errors,  may  require  increased  attentional  allocation,  and 
tends  to  contribute  most  to  the  task  difficulty. 


The  effect  of  increased  neuromuscular  lag  on  both  control  inputs 
also  appears  to  affect  the  lateral  axis  more  than  the  longitudinal  -  the 
azimuth  error  increasing  significantly  more  than  the  longitudinal.  As 
would  be  expected,  all  three  angular  rates  decreased,  but  these  were  higher 
than  the  simulated  values  for  all  three  models  -  thus  indicating  perhaps 
that  neuromuscular  lag  greater  than  0.2  seconds  is  warranted. 

TABLE  13 

SIMULATION  VS.  BASELINE  MODEL 


RMS  PERFORM 

ANCE 

Tracking  Error  (deg) 

£E1  eAz 

Lead  Angles  (deg) 

AE1  AAz 

Simulated 

Model 

1.72  2.58 

0.74  0.34 

Target  Velocities  (ft/sec) 
wT  vT 

Euler  Angles  (deg) 

e  <j> 

Simulated 

Model 

48.  34. 

37.  15. 

1.09  4.18 

0.74  2.29 

Velocity  Angles  (deg) 

a  6 

Angular  Rates  (deg/sec) 
p  q  r 

Simulated 

Model 

1.40  0.20 

1.03  0.43 

- 

4.99  1.72  0.92 

6.07  3.04  1.72 

Control  Inputs  (in) 

6r  <5p 

tST  mST  PED 

Simulated 

Model 

0.27  0.22  0.22 

0.28  0.23  _ 

TABLE  14 


THREE  BASIC  MODELS-RMS  PERFORMANCE 


Tracki ng 

rEl 

{ 

Error  (deg) 

i 

eAz 

Lead  Angles  (deg) 

*E1  AAz 

Basel i ne 

1.43 

0.80  j 

0.74 

0.34 

Incr.  Obs.  Err. 

1.60 

1.32 

0.97 

0.52 

Obs.  Err.  &  Lag 

1.66 

1.95 

0.34 

0.69 

Target  Velocities  (ft/sec)  ; 

Euler 

Angles  (deg) 

WT 

VT 

0 

<t> 

Basel i ne 

37. 

"t - 

15. 

0.74 

2.29 

Incr.  Obs.  Err. 

44. 

22. 

0.97 

3.32 

Obs.  Err.  &  Lag 

21. 

29. 

1.38 

4.70 

Velocity  Angles  (deg) 

Angular  Rates  (deg/sec) 

U 

a 

P 

q  r 

Basel i ne 

1.03 

0.43 

6.07  : 

1.  04  1.72 

Incr.  Obs.  Err. 

1.09 

0.47 

6.65  ' 

1.09  1.78 

Obs.  Err.  &  Lag 

1.15 

0.40 

6.30  \ 

1.58  1.38 

Control 

Inputs  (in) 

Motor  Noise  Ratio  (dB) 

* 

iS.  i 

est 

ast 

est 

ast 

Basel i ne 

0.28 

0.23 

-10 

-12.3 

Incr.  Obs.  Err. 

0.28 

0.27 

-10 

-13.6 

Obs.  Err.  &  Lag 

0.24 

0.25 

-  8 

-11.1 

Furthermore,  though  neuromuscular  lag  time  constants  on  the  order  of  0.1 
seconds  for  a  single-axis  task  and  a  force-sensitive,  side-arm  control 
stick  may  be  reasonable,  higher  time  constants  are  appropriate  for  a  movabl 
center-stick  control,  and  even  larger  time  constants,  approaching  one 
second,  for  rudder  pedals  in  this  multi -axis  task. 

Finally,  none  of  these  models  were  able  to  match  the  much  higher 
simulated  lead  angles  (a^  and  a^).  Furthermore,  none  of  the  other  models 
to  be  discussed  below  yielded  significantly  larger  rms  lead  angles.  The 
answer  probably  lies  in  the  fact  that  the  simulation  run,  although  almost 
representative  of  a  true  steady-state  aircraft  flight  condition,  did  not 
result  in  a  constant  range  to  the  target.  The  range  increased  from  the 
smallest  value  of  1800  feet,  occurring  early  in  the  simulation,  to  over 
3000  feet  at  the  end  of  the  40  second  run.  Consequently,  the  projectile 
time  of  flight  to  the  target  increased  proportionally.  Since  the  lead 
angle  tends  to  be  proportional  to  the  time  of  flight,  the  lead  angles  at 
the  larger  ranges  later  in  the  simulation  run  are  significantly  higher 
than  those  occurring  earlier  in  the  run,  which  are  closer  to  the  modeled 
constant  range  of  2000  feet.  Hence  the  rms  lead  angles  obtained  from  the 
simulation  would  be  expected  to  be  larger  than  those  prediced  via  the  model 

The  fact  that  the  largest  deviation  is  in  the  azimuth  lead  annle  is 
consistent  with  this  reasoning  as  well.  At  the  highly  banked  flight 
condition  considered  (<j>  -  75°),  the  largest  component  of  the  gravity 
vector  will  be  in  the  azimuthal  direction  rather  than  in  the  elevational. 
Hence,  the  azimuth  lead  angle  will  be  more  effected  by  the  "gravity  drop" 
of  the  projectile,  and  therefore  experience  the  largest  increase  as  the 
time  of  flight  increases. 

As  a  result  of  these  considerations,  we  shall  consider  the  aircraft 


states  as  the  most  important  for  comparison  purposes,  followed  by  the 
tracking  errors,  and  finally  the  lead  angles.  As  a  final  note  regarding 
the  matching  of  rms  performance,  the  simulation  run  represents  a  single 
case  involving  a  single  pilot,  hence  is  in  effect  a  single  data  point.  On 
the  other  hand,  the  statistical  results  from  the  model  represents  ensemble 
data  over  a  large  sample  Therefore,  too  much  emphasis  on  precise 

matching  is  foolish. 

With  these  initial  results  in  hand,  we  shall  now  address  the  issue 
of  the  inclusion  of  the  third  control  input,  the  rudder  pedal. 

Adding  the  third  control  parameter  to  the  model,  utilizing  a  pre¬ 
selected  neuromuscular  time  constant,  must  be  done  with  care.  The  two 
lateral-di rectional  controls  (aileron  and  rudder)  do  not  yield  pure  modal 
responses.  That  is,  an  aileron  deflection  produces  both  roll  and  yaw,  as 
is  the  case  with  the  rudder  as  well.  Recalling  the  control  equation  for 
the  pilot  model  as 

Vp  =  -V  -  u 

the  result  of  this  coupled  response  is  a  nondiagonal  T^  matrix.  Therefore, 
the  eigenvalues  of  T^  represent  the  effective  neuromuscular  lag  as  before, 
but  this  lag  is  for  a  combination  of  control  inputs  described  by  the 
eigenvectors  of  T^.  Therefore,  the  eigenvalues  of  T^  will  be  used  for  the 
lag  time  constants,  realizing  that  rather  than  representing  a  pure  aileron 
input,  for  example,  the  appropriate  time  constant  (or  eigenvalue)  represents 
primari ly  aileron  input  along  with  some  lesser  rudder  input. 

The  model  results  with  rudder  included  are  given  in  Table  15, 


along  with  the  simulated  results.  Note  that  the  three-control  model 
includes  the  higher  observation  noises  (fractional  attention  of  0.125  and 
angular  rate  thresholds  of  0.18  deg/sec)  as  well  as  the  larger  elevator 


TABLE  15 


BASIC  THREE-CONTROL 

MODEL  PERFORMANCE 

Tracking 

Error  (deg) 

Lead  Angles  (deg) 

£E1 

AE1 

AAz 

Simulated 

1.09 

0.97 

1.72 

2.58 

Model 

1.66 

1.03 

0.97 

0.47 

Target  Velocities  (ft/sec) 

Euler  Angles  (deg) 

WT 

VT 

e 

Simulated 

43. 

34. 

1.09 

4.18 

Model 

46. 

20.  ! 

0.92 

2.58 

Velocity  Angles  (deg) 

Angular 

Rates  (deg/sec) 

a 

6  i 

P 

q  r 

Simulated 

1.40 

0.20 

4.99 

1.72  0.92 

Model 

0.97 

0.63 

.  j 

5.44 

2.52  1.32 

Control 

Inputs  (in) 

Neuromotor  Lags  (sec) 

l 

5e 

tST 

5a  5r 

rtST  PED 

'NE 

XNA  tNR 

Simulated 

0.27 

0.22  0.22 

-- 

-- 

Model 

0.22 

0.39  0.28 

0.3 

0.3  1.1 
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lag  (0.3  seconds).  And  therefore  is  similar  to  the  third  model  for  which 
the  results  are  qiven  in  Table  14. 

First  note  that  much  higher  Ian  on  the  rudder  input  was  reauired. 
Attempting  to  use  a  lag  time  constant  of  around  0.3  seconds  for  all  three 
controls  resulted  in  lateral  directional  control  inputs  an  order  of 
magnitude  higher  than  those  shown.  Also,  one  would  expect  a  higher  laq 
for  pedal  controls  than  for  stick. 

Note  that  most  of  the  rms  values  for  this  three-control  model  are 
closer  to  the  simulated  results  than  the  similar  two-control  model  (the 
third  set  of  results  given  in  Table  14).  The  exceptions  are  in  the  much 
higher  side-slip  (t)  and  much  lower  bank  angle  (<j>).  We'll  return  to  this 
later. 

First  consider  the  results  depicted  in  Figure  12,  which  qives  the 
aileron-to-azimuth-error  pilot  describina  functions  obtained  from  the  two- 
and  three-control  models.  These  describing  functions  are  equivalent  to 
those  given  in  the  pitch-tracking  analysis  of  the  previous  section  in  that 
they  are  comparable  to  those  for  a  compensatory  task.  Note  that  the  low 
frequency  phase  of  the  two-control -model  result  is  zero,  while  that  of 
the  three  control  model  is  near  -270  degrees.  This  fact,  along  with  the 
comparison  of  the  low  frequency  slopes  of  the  magnitude  plots  (zero  for 
two  controls,  -20dB  per  decade  for  three  control  syndicates  that  the  de¬ 
scribing  function  qains  are  of  opposite  sign!  Which  one  is  correct? 

We  will  only  be  able  to  answer  part  of  that  question  here.  One  sign 
on  the  describing  function  corresponds  to  a  piloting  strateqy  of  rolling 
to  null  azimuth  error.  This  is  the  only  possibility  if  aileron  is  the 
only  control  input.  However,  with  both  control  inputs,  the  possibility  of 
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a  pointing  rather  than  a  rolling  strategy  is  possible,  with  "cross 
control"  inputs  being  used.  Inclusion  of  the  rudder  is  necessary,  then, 
to  allow  either  strategy  to  be  used,  and  the  existence  of  opposite 
sign  on  the  three-control  aileron  describing  function  would  represent 
the  effects  of  "pointing".  The  preliminary  conclusion  then,  would 
appear  that  inclusion  of  the  rudder  may  be  necessary  to  capture  the 
proper  piloting  strategy  "Which  strategy  is  correct?"  is  now  the  question. 

Now  let's  reconsider  the  rms  performance  of  the  three-control 
model  in  Table  15,  and  note  the  much  higher  value  for  rms  side-slip 
angle  (e)  and  the  lower  rms  bank  angle  ( <j> )  compared  to  the  previous  two- 
control  models  (Table  14).  The  fact  that  the  rms  azimuth  error  is  reduced 
in  the  three-control  case,  alonq  with  the  large  increase  in  side-slip 
and  lower  bank  angle  further  supports  the  presence  of  a  "pointinq"  rather 
than  "rolling"  strateqy.  Here  the  pointinq  may  be  accomplished  with 
rudder  input,  however.  But  the  fact  that  the  simulation  results  show 
much  lower  side  slip  indicate  that  the  rolling  strategy  was  actually  used. 
This  leads  to  the  conclusion  that  the  original  objective  function 
(weighting  only  displayed  errors  and  lead  angles)  is  insufficient  to 
correctly  model  the  pilot's  strategy  in  the  lateral -directional  axis. 

To  determine  the  required  modification  to  the  objective  function, 
we  must  consider  what  additions  to  the  pilot’s  observation  vector  are 
appropriate.  (Comments  from  a  test  pilot,  incidentally,  were  helpful  in 
this  regard.)  In  addition  to  lead  angles  and  tracking  error  (or  line  of 
sight),  the  pilot  also  sees  the  target  itself  through  the  heads-up  dis¬ 
play.  Consequently,  target  bank  angle  and  heading,  or  more  precisely, 
relative  bank  angle  and  heading  are  observed,  although  in  this  in-plane 
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maneuver,  the  relative  heading  was  considered  difficult  to  observe 
accurately.  However,  relative  bank  angle  was  observable,  and  in¬ 
tuitively,  if  nothing  else,  would  appear  to  be  important.  Furthermore, 
in  the  purely  single  axis  experiments  of  Harvey,  target  and  attacker 
bank  angle  were  both  always  zero,  hence  no  importance  of  bank  angle 
could  be  ascertained  from  his  analysis. 

This  raises  another  key  point  not  previously  considered.  Harvey's 
experiments  were  fixed  base,  while  in  the  Air  Force  LAMARS  facility, 

motion  cues,  and  particularly  accelerations,  were  most  definately  available. 

n  21 

Furthermore,  clear  evidence  is  available1  J  that  pilots  avoid  lateral 

accelerations,  while  vertical  accelerations  of  reasonable  magnitude  (e.q. 

less  than  5  to  6  g's)  are  not  as  highly  disliked,  or  at  least  pilots  are 

more  conditioned  to  the  latter  type.  Consequently,  inclusion  of  relative 

bank  angle  and  lateral  acceleration  in  the  observation  vector  is  clearly 

reasonable,  and  weighting  both  of  these  parameters  in  the  objective  function 

is  consistent  with  this  entire  discussion  of  piloting  strategy. 

To  complete  the  model  then  the  observation  threshold  on  relative 

bank  angle  was  taken  to  be  5  degrees,  based  on  the  pilot's  perception 

thresholds  of  0.05  degrees  at  the  pilot's  eye,  tracking  ranqe  of  2000  ft., 

and  target  wing  span  of  approximate! y  40  feet.  The  perception  threshold 

on  lateral  acceleration  was  assumed  to  be  0.4  ft/second  or  0.012  g's, 

f  1 3"1 

roughly  consistent  with  some  earlier  studies:  J  Finally,  using  an 
analysis  similar  to  that  above  for  bank  anqle,  it  was  determined  that  the 
threshold  on  relative  bank  angle  rate  is  approximately  10  deg/sec.  But 
since  this  is  much  larger  than  the  rms  roll  rates  (p)  encountered,  no  useful 
information  could  be  gained  by  including  roll  rate  in  the  observation  vector. 


Inclusion  of  the  above  parameters  results  in  the  multi-axis  pilot 
model  summarized  in  Table  16  .  Note  that  since  the  target  has  a  constant 
bank  angle,  observing  and  weighting  relative  bank  angle  is  equivalent  to 
observing  and  weighting  the  attacker’s  bank  angle  in  this  case.  Finally, 
the  attentional  allocation  has  been  adjusted  to  further  improve  the  model 
results,  and  these  allocations  reflect  the  relative  difficulty,  or  work¬ 
load,  between  the  two  axes. 

With  the  above  model,  the  rms  performance  shown  in  Table  17  was 
obtained.  The  parameters  marked  with  the  double  asterisk  indicate  this 
value  was  the  closest  to  the  simulated  value  of  all  models  considered, 
while  the  single  asterisk  denotes  one  other  model  considered  resulted  in  a 
closer  match  for  that  parameter.  Considering  all  the  factors  discussed, 
it  is  felt  that  this  model  is  quite  justified. 


SYSTEM  AUGMENTATION 

With  this  final  pilot  model  established,  we  may  proceed  with  the  system 
augmentation.  Reiterating  the  methodology,  with  augmentation,  the  system 
dynamics  are  represented  as 


Ax 


+  BPUP  +  BAuA 


+  w 


where  Up  and  ua(jg  are  the  vectors  of  control  inputs  of  the  pilot  and  aug¬ 
mentation,  respecti vely,  while  w  represents  any  system  disturbances.  The 
control  u_M„  is  chosen,  then,  to  minimize  the  objective  function. 

a  Uy 


aug 


=  +  E 


lim  i 

T-**=  1 


rT 


Fu 

aug  aug 


u' 


dt 
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SB 


Observation  Vector,  y  -  C E ^ -j »  £E1*  cAz’  eAz*  ^E1  *  ^E1  ’  ^Az’  ^Az*  S’  ^rel^ 
Objective  Function  Weights,  Q  =  16. ,  Q*  =  1 . ,  Q?  =  4. ,  Q  =  .007,  Q  =  8. 


Observation  Thresholds,  T  =  T  *  .05  deg,  T*  =  T*  =  .18  deg/sec 

T  =  .4  ft/sec^,  T  =  5  deg 
ay  * 


Fractional  Attention  Allocations,  f.  = 

.05  on 

fi  = 

.15  on  eAz,  eAz,  xAz,  xAz 

fi  = 

.1  on  a  and  <*> 

y 

Observation  Noise  Ratio  -20dB  for  full  attention 


Observation  Delay,  t  =  0.2  sec 


Neuromuscular  Lag,  tn  =  .33  on  6  (obj.  func.  wt.  =  .05) 

n  bst 

tn  =  .23  on  6.  (obj.  func.  wt.  =  .10) 
st 

=  .62  on  6r  (obj.  func.  wt.  =  .12) 
ped 


2 

Motor  Noise  Ratios,  Vp  =  ^pu^u 

p  =  .05  (-8  dB)  on  5F 

u  St 

p  =  .02  (-12  dB)  on  6. 
u  st 

p  =  .05  (-8  dB)  on  6D 
u  ped 


TABLE  17 

FINAL  MODEL  PERFORMANCE 


Tracking  Error  (deg) 


Lead  Angles  (deg) 


Simulated  1.09  0.97 

Final  Model  1.78  1.72 


Target  Velocities  (ft/sec) 


Wj  v. 


48. 

34. 

63. 

25. 

Velocity  Angles  (deg) 

a  3 


Simulated  1.40  0.20 

Final  Model  0.97  0.15** 


Control  Inputs  (in) 

SE  ^A  5R 
st  st  ped 

Simulated  0.27 
Final  Model  0.23 


Lateral  Acceleration 
(ft/sec2) 


1.72 

2.58 

1 .38** 

0.63* 

Euler  Angles 

(deg) 

e 

♦ 

1.09 

4.18 

0.97* 

3.55* 

Angular  Rates 

(deg/sec) 

p  q 

r 

4.99  1.72 

0.92 

5.21**  2.58** 

0.57** 

Neuromotor  Lags  (sec) 

tne  tna  tnr 

.33  .28  .62 

Neuromotor  Noise  Ratios 
(dB) 


Simulated 
Final  Model 


■ 

.6 

.7; 

-12  -8 
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where  Jp  is  the  pilot's  objective  function  as  determined  previously  in 
the  pilot/vehicle  analysis.  The  optimal  control  law  is  then  (as  shown 
previously) 

uaug  =  *  8IKS2Up) 

where  Kfl  and  Kfl  are  determined  by  solving  the  pilot  and  augmentation 
12 

Riccati  equations  simultaneously.  The  augmented  piloted  system  may  then 
be  represented 

*  *  <A  -  *  ” 

or  (A  -  B^F  )  =  (A  -  B^Kx)  =  augmented  plant  matrix 

(Bp  -  B^F  )  =  (Bp  -  B^K^)  =  augmented  control  matrix 

(e.g.,  stick  gains  and  cross  feeds) 

The  pilot's  control  law  remains 
Up  =  -G-VpKp  x  -  G” 1  BpKp  Up 

except  the  Riccati  gains  KD  and  KD  now  reflect  the  presence  of  augmentation 

1  2 

as  they  are  solved  simultaneously  with  the  augmentation  problem  using  the 
augmented  plant  matrices. 

The  open-  and  closed -loop  root  loci  as  augmentation  level  increases 
from  unaugmented  (denoted  by  U )  to  higher  levels  of  augmentation  (denoted 
1-4)  are  shown  in  Figures  13  and  14  .  The  first  figure  shows  the  (primarily) 

longitudinal  eigenvalues,  recognizing  this  system  actually  has  coupled 
longitudinal  and  lateral -directional  modes,  with  the  second  figure  depicting 
the  (primarily)  lateral-directional  roots.  As  with  the  pitch  tracking 
results  of  the  previous  section,  the  augmented  system  roots  are  denoted  as 
"X"  and  the  piloted  closed-loop  roots  are  denoted  as  "0". 


Examination  of  the  longitudinal  root  locus  leads  to  the  conclusion 


that  the  trends  are  changed  only  slightly  from  those  obtained  in  the 
single-axis  pitch-tracking  analysis.  The  two  low  frequency  roots  from 
the  geometry  and  kinematics  of  the  problem  were  at  the  origin  for  the  wings- 
level  pitch  analysis,  but  still  move  quickly  toward  the  invariant,  low- 
frequency  pair  of  piloted  closed-loop  roots  with  augmentation.  Recall 
this  probably  indicates  the  exi stance  of  a  pair  of  pseudo-zeros  or 
transmission  zeros  at  this  point.  The  trend  in  the  vehicle  short-period 
roots  is  also  similar  to  previous  results,  with  the  exception  that  short 
period  damping  is  increased  much  more  due  to  augmentation.  Finally,  the 
proximity  of  the  open-loop  roots  to  the  closed-loop  roots  for  the  unaug¬ 
mented  plant  ( U)  and  level-two  augmented  plant  (2)  indicate  less  compensation 
being  introduced  by  the  pilot,  primarily  due  to  his  reduced  bandwidth 
(xn  =  .3  instead  of  xn  =  .1)  in  this  more  complex  task. 

The  root  loci  for  the  lateral-directional  axis  will  now  be  considered. 
The  unaugmented  system  has  two  eigenvalues  near  +  .2j  associated  with 
the  azimuth  portion  of  the  geometry  and  kinematics  of  the  engagment,  the 
spiral  and  roll  subsidence  vehicle  roots  at  the  origin  and  near  -1  per 
second,  respectively,  and  the  ditch  roll  roots  at  -.4  +  4.7j  rad/sec. 

With  augmentation  now,  the  damping  and  frequency  of  the  dutch  roll  are 
increased,  and  a  complex  roll/spiral  mode  is  formed  from  the  spiral  and 
roll  subsidence  roots.  The  two  low-frequency  roots  on  the  jw  axis  ultimately 
attach  to  the  real  axis  near  -.6  and  become  two  real  roots.  The  first  of 
these  approaches  the  invariant  closed-loop  root  at  -.5/sec,  the  other 
moves  to  the  left  on  the  real  axis.  This  latter  root,  incidently,  appears 
to  be  attracted  to  the  pilot  closed-loop  roots  near  -1/sec  which  also 
approach  the  real  axis  with  increased  augmentation.  Finally,  associated 
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with  the  dutch-roll  and  the  roll-spiral  are  two  pairs  of  piloted  closed- 
loop  roots  very  close  to  their  associated  unagumented  branches. 

Except  for  the  roll-spiral  mode  without  augmentation,  the  pilot  appears 
to  introduce  relatively  little  compensation  in  terms  of  eigenvalue 
repositioning,  again  due  to  the  complexity  of  this  multi-variable  task 
and  his  limited  bandwidth.  The  relatively  large  importance  of  this  roll/ 
spiral  mode  could  be  inferred  from  this  result. 

The  level-two  augmentation  was  selected  for  evaluation  and  it  appears  to 
represent  a  medium -to  high-level  of  augmentation.  (The  gains  and  Ku  are 
given  in  Table  18  for  all  augmentation  levels).  The  frequency  response  for 
the  system  and  pilot  describing  functions  are  presented  in  Figures  15  to  20 
for  unaugmented  and  augmented  systems. 

The  longitudinal  results  shown  in  Figure  15  indicate  the  near  cancellation 
of  the  low  frequency  eigenvalues  due  to  closure  of  the  tracking  loop,  and  in¬ 
creased  damping  of  the  short-period  mode  (near  3.5  rad/sec)  is  evident.  The 
pilot's  describing  function  shown  in  Figure  16  shows  little  change  due  to  aug¬ 
mentation. 


The  aileron  frequency  responses  are  given  in  Figures  17  and  18  .  The  phase 
results  for  the  unaugmented  vehicle  clearly  show  the  presence  of  the  spiral 
root  at  the  origin,  as  well  as  the  two  low-frequency  kinematic  roots.  These 
three  roots  appear  to  be  effectively  cancelled  with  augmentation. 

The  pilot's  aileron  describing  function  indicates  almost  a  constant  aileron- 


to-error-rate  tracking  characteristic,  or 


yp(s)  =  V5) 


With  augmentation,  the  slope  of  the  magnitude  curve  is  reduced,  implying  lesser 


emphasis  on  rate.  This  should  indicate  less  pilot  workload  as  measured  by  lead 


requirements. 
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TABLE  18 


-0.001  0.056  -0.002  0.040  -0.002  0.017 


AUGMENTED 


Finally,  the  rudder  frequency  responses  are  shown  in  Figures 
19  and  20.  As  with  the  aileron  results,  the  pilot's  rudder  describing 
function  indicates  that  the  pilot  is  providing  turn  coordination  and 
regulating  azimuth  error  with  rudder  input. 

To  see  the  effects  of  the  augmentation  on  the  statistical  perfor¬ 
mance,  consider  the  results  in  Table  19.  As  anticipated,  the  rms  tracking 
errors  are  reduced,  and  it  is  noted  that  the  azimuth  error  has  been  re¬ 
duced  more  than  the  elevation  error.  Since  the  lateral  directional 
axis  is  considered  the  more  difficult,  it  appears  that  the  augmentation 
has  significantly  improved  the  tracking  performance  as  predicted  by 
linear  analysis. 

Note,  finally,  that  the  variance  of  the  additive  motor  noise  was  the 
same  for  the  augmented  system  analysis  as  for  unaugmented.  This  results 
in  a  higher  noise- to-signal  ratio  for  the  motor  noise  in  the  augmented 
case,  and  should  yield  conservative  results  for  the  performance  pre¬ 
dictions.  In  addition,  the  piloted  predictions  were  made  assuminq  the 
pilot  used  all  three  controls.  However,  the  augmentation  gains  were  de¬ 
termined  using  the  assumption  that  no  pilot  rudder  input  would  be  used 
to  reduce  pilot  workload  (hence  no  gains  on  6^  in  Table  18).  If  no  pilot 
rudder  was  actually  evaluated,  rms  results  would  be  further  improved. 
Simulations  are  required  in  any  future  efforts  to  substantiate  these  re¬ 
sults. 
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TABLE  19 

AUGMENTED  SYSTEM 
RMS  PERFORMANCE 


Tracking  Errror  (deg) 

Lead  Angles  (deg) 

eEl  eAz 

XE1 

XAz 

Unaugmented 

1.78  1.72 

1.38 

0.63 

Augmented 

0.54  0.20 

0.10 

Target  Velocities  (ft/sec) 

Euler  Angles  (deg) 

WT  VT 

e 

$ 

Unaugmented 

63.  25. 

0.97 

3.55 

Augmented 

15.  4. 

0.17 

0.48 

Velocity  Angles  (deg) 

Angular 

Rates  (deg/sec) 

a  6 

P 

q  r 

Unaugmented 

0.97  0.15 

5.21 

2.58  0.57 

Augmented 

0.36  0.06 

1 .20 

0.97  0.21 

Control  Inputs  (in) 

Neuromotor  Lags  (sec) 

5e  5a  5r 

tst  rtst  ped 

Tna  \ 

Unaugmented 

0.23  0.15  0.10 

0.33 

0.28  0.62 

Augmented 

0.19  0.11  0.06 

0.33 

0.28  0.62 

Lateral  Accelerational 

(ft/sec2) 

Neuromotor  Noise  Ratios 
(dB) 

ay 

N. 

6E 

N,  N 

6a  6r 

Unaugmented 

0.78 

-8. 

-12.  -8. 

Augmented 

0.25 

-7. 

-8.  -4. 
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APPENDIX  A 


EQUATIONS  FOR  THE  RELATIVE  LINE  OF  SIGHT 

Vector  Equation  Development 

Clearly,  the  angle(s)  to  the  line  of  sight  (LOS)  are  required  for  our 
analysis,  and  as  usual,  linearized  perturbation  equations  are  ultimately 
required.  The  equations  to  be  developed  in  this  section  will  be  derived 
for  the  angles,  as  shown  in  the  figure  below,  in  the  vehicles  gun  coordi¬ 
nate  system,  X  Y  Z  . 

9  9  0 


Now  the  kinematics  of  the  line  of  sight  unit  vector,  s,  are  described  by 
’target  '  Stacker  ’  5T<0i>  =  +  « 

where  D  is  the  scalar  distance  from  the  attacker  to  the  target. 


But  the  time  rate  of  change  of  a  unit  vector  is  the  cross  product  of  that 
vector  with  its  rate  of  rotation  vector  in  inertial  space,  '—■>  or 

P 

?  =  ?  X  s 

So  6  is  the  rate  of  rotation  of  the  LOS  in  the  inertial  (e.g.  earth- 
fixed)  coordinate  system. 
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Taking  the  cross-product  of  Equation  A1  with  s  yields 

•  •  • 

?  X  (VT  -  VA)  =  s  X  Ds  +  s  X  (B  x  s)  D  =  Ds  X  (f  X  s) 

And  from  the  vector  triple  product  identity 

•  •  • 

s  X  (B  X  s)  =  B  -  (B  •  s)  s 

Now,  instead  of  the  inertial  LOS  rate  8,  we  desire  the  LOS  rate  in 
the  attacker's  gun  coordinate  system  or  b"  ,  which  is 

gun  p  g 

where  in  the  rotation  rate  vector  of  the  gun  or  aircraft  system.  Sub¬ 
stituting  into  the  above  yields 


V  *  "g  ■  [<V  +  V  '  s>  ■  B  5  x  <vT  -  ».) 


or 


=  ^  s  K  <VT  '  V  '  “q  +  U  •  s)s  +  (u>  •  s)s 


gun 


Mow  the  term  ( 6gun  •  s’)  =  0  since  the  rotation  vector  and  the  line  of 
sight  in  the  gun  coordinate  are  orthogonal. 

The  necessary  vector  equation  is  then 


8  -  i[s  X(VT  -  VJ]  -  uj  +  (w  •  s)s 

gun  DL  '  T  A  J  g  '  g 


(A2 ) 


From  the  previous  figure,  we  can  define  the  vectors  as  follows,  using 
small  angle  approximations 

?  =  V  BAzJg  +  eEl¥g 

®gun  “  '6E1JK  +  r:Azkg 

(D  =Pi  +  Q  j  +RK 
g  g  g  MgJg  g  g 

V.  =  u.T  +  VJ  +  w.k 
A  A  g  AJg  A  g 

V,  =  UT  T  +  VT  T  +  WT  ¥ 

I  I  q  I  0  T  a 

g  y  g  y  g  y 
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(It's  important  to  note  that  Vy  and  the  components  here  defined  are  the 
velocity  components  of  the  target  in  the  attacker's  gun  coordinate  system.) 
Carrying  out  the  operations  in  Equation  A2,  and  assuming  w  •  s  ;  P^,  we 
have  the  desired  relations 

»E1  ■  n  [<“xg  -  V  -  eE1(UTg  -  V]  +  Qg  ■  Yaz 

»AZ  *  F  t'VTg  -  V  -  »Az<UTg  "  Us>J  -  R,  +  YeI 

To  obtain  the  equation  for  the  rate  of  change  of  range,  or  D,  we  can 

simply  take  the  dot  product  of  Equation  A1  with  s  to  obtain 

•  • 

(VT  -  VA)  •  ?  =  D(s  •  s)  +  D(?  X  s)  •  s 
so 

D  =  (VT  -  VA)  *  ? 

D  =  (UT  -  U  )  +  (VT  -  V  )ba  ♦  (WT  -  W  )b 
g  g  g 


“i 


Seal  or  Equations  in  Stability  Axes 

Since  all  our  equations  are  ultimately  expressed  in  the  vehicle  stability 
axes,  we  need  to  express  the  above  equations  in  terms  of  stability-axes 
variables  rather  than  gun-coordinate  variables. 

Assuming  the  gun  axis  to  be  aligned  with  the  vehicle  body  reference 
axis,  we  havf  the  following  relations  between  gun  ( • )  and  stability  (*)s 
variables 

'lg  =  Us  cos  a-j  -  W$  sin  a-j 


Wg  =  Ws  cos  +  Us  sin 
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Pg  =  Ps  cos  a-|  -  R$  sin  a-j 


Q  =  Q 

wg  vs 

Rg  =  Rg  cos  a-j  +  P$  sin  aj 


The  resulting  equations  are 


D  =  [(UT  -  U  )  cos  -  (WT 
s  s 

*  SAz<VTs  -  V 


W$'  sin  a-j ] 


+  6E1[(WT  -  Ws)  cos  +  (UT  -  U$)  sin  c^] 


6E1  =  D^WT  "  Ws^  C0S  “l  +  ^UT  "  V  Sin  al 
s  s  ' 

-  BE-j{(UT  -  Us)  cos  a-j  -  (WT  -  W$)  sin  a1 } ] 
s  s  ’ 

+  Qs  -  6Az(ps  cos  -  Rs  Sin  Hj) 

Us)  cos  -  (WT  -  W$)sin  a-j}] 
(R$  cos  a-j  +  Ps  sin  a-j)  +  eEL(Ps  cos  a-j  -  R^  sin  r;-j  ) 


(A3) 


BAz  "  D^VT  *  V  ‘  eAz{(UT 

s  s 


where  U-j-  ,  V,  ,  W_  are  now  the  components  of  the  target  velocity  in  the 
s  s  s 

attacker's  stability  axis  and  a-j  is  the  attacker's  trim  angle  of  attack 
(or  the  trim  angle  between  the  gun  axis  and  stability  axis). 


The  Perturbation  Equation 

Now  to  linearize  the  above  equations,  we  will  define  the  motion  variables 
in  terms  of  a  constant,  steady-state  quantity  (O-j  and  a  time-varying  per¬ 
turbation  about  the  steady-state  variable.  So,  let 

us  5  U1  +  u  ps  =  P1  +  P 

W  v  Qs  »  Q,  *  C 

W$  =  W-j  +  w  R$  =  R^  +  r 


D  =  Dj  +  d 
BEL  =  BELj  +  BE1 

8Az  =  BAz]  +  BA 
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Inserting  these  into  Equations  A3  will  yield  two  sets  of  equations, 
one  for  the  steady  state  values  and  one  for  the  perturbation  quantities. 

For  example,  taking  the  equation  we  have 

(D1  +  d)(ftE1  +  eE1)  =  [(WT  -  W] )  cos  +  (U^  -  U]  )sin  «j] 

+[(wT  -  w)  COS  a1  +  (uT  -  u)  sin  a-j  ] 

-(bEi  +  eE1)[(U-.  -  U] )  cos  a}  -  (WT^  -  W-j  )sin  (A4) 

+( u j  -  u)  cos  a-j  -  (Wj  -  w)  sin 
+  (D-j  +  d)(Q1  +  q) 

-  (D-j  +  d ) ( Bftz ^  +  bAz)[(P1  +  p)cos  -  (Rj  +  r)sin  Oj] 


Now  the  steady  state  (•).  values  certainly  must  obey  the  original  equation, 
so 


o  =  [(w, 

'l 


+  D 


1^1 


-  W-j  )  COS  a-j  +  (Uj. 
C(UT  -  Uj)  cos  aj  - 

D1  BAZjLPlCOS  al 


-  Uj)  sin  aj] 

(WT  -  Wj )  sin  <Xj] 

-  Rj  sin  «j] 


We  may,  as  a  result,  subtract  this  relation  out  of  the  complete,  original 
relation  (Eqn  A4).  In  addition.,  under  the  assumption  that  the  perturbation 
variables  are  small,  we  will  drop  higher  order  terms  (e.g.,  products  of  per¬ 
turbation  quantities).  Finally,  noting  that  BEE,  =  0  by  definition,  we  have 

DleEl  =  ^WT  "  cos  al  +  (UT  ”  u)  sin  al^ 

-  Bf1  [(u,  -  u)  cos  a,  -  (wT  -  w)  sin  a,] 

tll  ‘l  1 

-  6E1[(UT  -  Uj)  cos  0j  -  (WT  -  Wj)  sin  Oj]  (cont'd) 
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+  D^q  +  Q^d  -  d  B^z  (P^  cos  a ^  -  R-j  sin  ) 

-  BAz  D1(P1  cos  a-j  -  R-|  sin  (p  cos  a-j-r  sina.j 

This  equation  is  now  linear  in  terms  of  the  perturbation  variables, 
and  the  steady-state  variables  are  constant  by  definition. 

In  like  manner,  the  remaining  two  equations  in  A3  result  in 

D-j  =  0  =[(UT  -  U})  cos  0]  -  (W^  -  sin  c^]  + 

+  BE1  C(Wt  -  W1 )  cos  a-|  +  (Uj  -  U-j)  sin  o^] 

d  =  [(uT  -  u)  cos  o1  -  (wy  -  w)  sin  o^]  +  BAz(Vj  -  V-| ) 

+  BAz^(vt  -  v)  +  eE-,[(Wy  -  W] )  cos  a1  +  (Uy^-U^sin  a]  ] 

+  Be1  [(Wy  -  w)  COS  a-j  +  (Uy  -  u)  sin  ] 

BA2,"  °’b4<vt,  •  »,)  -  BA7«UT,-  V  cos  «1  -  (“t,  -  VSt"  V1 

-(R^  cos  a-j  +  P1  sin  (P^  cos  -  R1  sin  ) 

°l^Az  =  ^VT  "  “  hz^uJ  "  cos  al  ‘  ^WT  '  sin  al-^ 

"6Az*-(UT1  “  V  C0S  al  "  (WT  ‘  V  Sin  al^ 

-D^  (r  cos  a1  +  p  sin  a-j )  -  d  (R^  cos  +  P1  sin  ) 

+Dy BE1  (p  cos  a-j  -  r  sin  «i)  +  D1BE1  (P^  cos  a-j  -  R-j  sin  a-j ) 

+  d  B^i  (P.j  cos  a-j  -  R.|  sin  a-| ) 


and 
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APPENDIX  B 


EQUATIONS  FOR  THE  SIGHT  LEAD  ANGLES 


Perfect  Director 

Again,  linearized  equations  are  needed  for  the  lead  angles  computed 
by  the  sights.  From  the  reference  ,  the  fundamental  equation  for  the 
lead  angle, I,  in  general,  is 

Vfl  -  ?  x  (VT  -  VA)  ♦  JvVA?gun  ♦  (?  X  Vlf  <B1> 


where,  referring  to  the  figure  and  assuming  the  sight  (pipper)  is  on  the 
target,  we  have 

?!V  XAzJg  +  XE1  *9 

or 


i  +  B.  j  +  8,~iK 
g  AzJg  El  g 

s  X  ’gun  _  XEl^gun  ~  XAzkgun 

average  projectile  velocity  in  still  air 

D^/T^.  -  V^,  D^.  =  distance  of  flight) 

projectile  time  of  flight 

attacker's  velocity 

V  -  V 

Jump  correction  factor,  =  _M _ f,  VM  is  projectile  muz¬ 


zle  velocity 


VA  +  VM 


’gun  X  *  =  " °^gun  +  Bkgun’  °  and  6  here  are  angle 
|VA' 

of  attack  and  sideslip 
target's  acceleration 


★ 

Hohwiesner,  W.,  Capt.  "Principles  of  Airborne  Fire  Control,"  USAF  Academy, 
Dept,  of  Astro,  and  Comp.  Sci.,  Dec.  1975. 
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To  obtain  the  scalar  component  equations,  define  the  following  vectors 


in  the  gun  coordinate  system 


where  as  in  the  LOS  equation  development,  it's  important  to  remember  that  the 
target's  velocity  and  acceleration  components  are  in  the  attacker's  gun 
coordinate  system.  Carrying  out  the  vector  operation  yields 


VeI  ■  [8E1(UT  -  V  -  <“t  -  “g»  -  JvV 
9  9 


+  (0nAT  ”  at  ^  f 
El  T  Tz  T 

g  g 


-Va2  -  -  V  -  eAZ<UTg  -  Ug)]  +  JVVA6 


+  (AIy  -  bAz\ 

g  g 


These  equations  may  be  used  to  calculate  lead  angles  assuming  the 
actual  1 ine-of-sight  angles  (e's)  and  target  acceleration  are  known.  This 
ideal  case  is  referred  to  as  a  perfect  director  sight. 

Now,  these  equations  will  be  expressed  in  terms  of  stability  axis 
velocities,  etc.  rather  than  gun  axis.  We  have  the  relations 
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U  =  U  COS  a,  -  W,  sin  a, 
g  s  is  1 

Vvs 


W  =  W  cos  a,  +  U  sin  a. 
g  s  is  1 


Ay  =  Ay  cos  <*i  -  Ay  sin 
Xg  Xs  Zs 


Ay  "  Ay 

Yg  Ys 


Ax  =  At  cos  a,  +  A  sin  a, 

Tz„  \  T, 

g  s  x$ 


Finally,  using  the  perturbation  technique  employed  in  the  previous 


appendix  we  introduce  the  relation 


XE1  =  LE1  +  XE1 


a  =  +  a 


AT  =  AT  +  3t 
TXs  ,X1  \ 


Introducing  these  perturbations  and  assuming  that  Vf,  Tf,  and  Jy  are 
constant  for  this  analysis  we  have  the  following  steady  state  relations 

VfLEl  =  BEi  C ( Uj  -  Uj )  COS  a1  -  (Wy  -  Sin  c^] 

- [ (Wy  -  W,)  COS  a-j  +  (Uy  -  U,)  sin  a.j] 

“JvVA,al  "  C0S  al  +  ATv  Sin  al 

1  6  zi  *i 

~BE1  (Ay  cos  “l  "  AT  sin  “l^ 

X1  Z1 


^fLAz  =  "  V  BAz^UT1  '  V  C0S  al  ”  ^WT1  "  V  S1n  “1^ 

+  JvVl  +  -  BAz(Ay  cos  a,  -  A  sin  a,] 

£  Tl  Al  L 1 


\ 
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Furthermore  the  perturbation  equations  are 


/fA^  =  Be1L(Uj  -  u)  cos  <*1  -  (w-j.  -  w)  sin  ctj] 

+  0^[(UT^  -  U1 )  cos  -  (Wy  -  W])  sin  a-| ] 

-  [(w-j-  -  w)  cos  a-j  +  ( u T  -  u)  sin  c^] 

-Tf(ax  cos  a,  +  aT  sin  a,) 

T  Tz  1  Tx  1 


-  SE-j(a-j-  cos  a1  -  aT  sin  a-, ) 


1  T. 


and 


-  8^  (Ay  COS  -  Ay  sin  a-j) 

X1  Z1 

-  jvV' 


'VfXAz  =  ^ VT  "  "  BAz^uT  '  cos  al  "  ^WT  "  w)  sin  ai^ 

"6Az^UT1  '  V  C0S  al  '  ^WT1  "  V  Sin 

+  Tf[aT  -  B.  (aT  cos  a,  -  aT  sin  a,) 

2  1 Y  x  1  z 

"BAz^AT  C0S  al  "  AT  Sin  al^ 


+  JyVAB" 


LCOS  Equation 

Returning  to  equation  B1  in  the  previous  section,  recall  that  the  per¬ 
fect  director  assumes  perfect  knowledge  of  the  1 ine-of-sight  rate  arid  target 
acceleration.  The  lead  computing  optical  sight  (LCOS)  system,  however,  uses 
attacker  rotation  rates  and  acceleration  to  approximate  the  above  parameters 


In  terms  of  the  line  of  sight  rate  in  the  attacker's  gun  coordinate 
system,  s’  we  have  from  the  cited  reference, as  well  as  the  previous 
appendix, 


but,  the x vector  defined  here  is  used  as  an  approximation  for  -Bgun  f°r  the 
LCOS,  and  where  the  ?  vector  is  defined  in  the  previous  appendix. 


In  this  case  then, 


and 


so 


Now  since  as  shown  in  previous  appendix 
s  X  (VT  -  VA)  =  D[i  -  {? •  s)  s] 

?  -  (?•?)?=  (ug  -  I)  -  [(ug  -  I)  -  Sis 

•  • 

s  X  (VT-  VA)  =  (u  -  T)-(wg  •  s)s  +  (I  .  s)s 


In  addition,  the  attacker  acceleration  Aft  is  substituted  in  place  of 
the  target's,  resulting  in  the  vector  equation  for  the  LCOS  (with  x  •  s 
equal  zero) 

V  =  D[(“g  "  '  (Vs)*3  +  JvVgun  +  (s  X  V  -? 


so,  with 


‘‘V  JA2Jg  +  JELkg 


u_  =  +  +  R„k„ 

g  g  g  g  g  g  g 


and 


S.  =  A  T  +  A  J  +  A  7 
A  x  g  ydg  z  g 
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wr 

■ 


we  have^the  following  two  scalar  differential  equations 


DXEL  "  -VfXEL  +  DQg  '  °(Pg  +  XAzQg  +  XELRg^XAz 


-JvV  +  ^XE1Ax  -  Az} 


-DX.  =  VrA.  +  D  R  -  D(P  +  A.  Q  +  An  R  )  ACI 
Az  f  Az  g  '  g  Az  Hg  EL  g'  EL 

+J  V.B  +  Tf(A  -  A.  A  ) 
v  A  y  Az  x' 


The  steady  state  equations  are  then,  using  the  usual  perturbation 


technique. 


Di4l  ■  0  ■  'Vaz  +  -  JvYl 

+  ^jf  [LCI  (A  cos  a.  -  A  sin  a.)  -(A  cos  a,+A  sin  a-,)] 

2  tL  X  *j  I  Zj  I  Z -j  I  X-j  I 

Di L^z ^ ( P-|  ai  -  ai )  L^zQ-|  ^  ^"El_  cos  a-j  + 


P]  sin  o^)] 


-D]LAz  =  0  =  VfLAz  +  D1(R1  cos  a]  +  P  sin  a] ) 


*  *  4[\,  -  LAz  <Ax,  cos  al  -  Az,  s1n  al>] 

_D1 Lel[ ( P1  cos  al  "  R1  sin  al  >  *  LAzQ1  *  LEL(R1  cos  “1  * 


P]  sin  ot1 )] 


Finally,  the  perturbation  equations  become 


D]  A^|^  -  "VfA^i  +  D^q  +  Q^d  -  OyVAa 


+  __f [LE-^(ax  cos  a-|  -  az  sin  ) 

+  A'  (A  cos  a.  -  A  sin  a.)  (cont’d) 

L.  I  X-j  I  Z-j  I 


-  (az  cos  dj  +  ax  sin  ctj )] 

-  dLAZMl  -  VazM1 

■  -  D1LazC(p  COS  -  r  sin  ^ )  +  X^Q,  +  L^q 
+  X^(R-j  cos  a-j  +  P-|  sin  ) 

+  l_EE(r  cos  ai  +  P  sin  al ^ 

‘Vaz  =  VfxAz  +  Vr  C0s  al  +  P  sin  V 

+  d(R1  cos  a-j  +  P1  Sin  c^)  +  JyVAB' 

+  Ift*y  -  lkz\  C°S  “1  '  *2,  Sl"  “l' 

-  La  (ax  cos  o-j  -  az  sin  c^)] 

-  dL£LM1  -  D1XE1M1 

-  D1LEL[(p  cos  a,  -  r  sin  a-, )  +  LAZq  +  X^Q, 

+  X'^  cos  0]  +  P1  sin  a})  +  LEL(r  cos  a,  +  P  sin  c^)] 


where  in  the  above  relations 


M-J  =  (P1  cos  a1  -  R-,  sin  )  +  LAzQ1 
+  Lel(R-j  cos  a 1  +  P1  sin  c^) 


APPENDIX  C 


TARGET  KINEMATICS 


In  the  foregoing  appendices,  the  equations  for  the  I i ne-of-sight 
rates  and  sight  lead  angles  were  developed  in  terms  of  the  attacker's 
velocities  and  of  the  target's  velocities  expressed  in  the  attacker 
coordinate  system.  In  this  appendix,  the  relations  governing  these  tar¬ 
get  kinematics  will  be  developed. 

Consider  the  target's  velocity  vector  expressed  in  the  attacker's 
gun  coordinates  as 


^T  ^T  nA  +  ^T^A 


"A 


Now,  the  time  rate  of  change  of  this  vector  with  respect  to  the  attacker's 
coordinates  is  simply 


=  UT  TA  +  VT  h  +  WT 


But  we  may  also  express  this  acceleration  vector  in  terms  of  the  time 
rate  of  change  of  the  target's  velocity  vector  relative  to  the  target's  own 
coordinate  system,  plus  the  effect  of  the  relative  motion  between  the  two 
coordinate  systems. 


or 


T  +  urel 


(Cl) 


where  V 


Here,  then,  u>  .j  is  the  rate  of  rotation  of  the  target's  coordinate 
system  relative  to  the  attacker's,  or  in  terms  of  the  two  aircraft  angular 
rates  relative  to  inertial  coordinates 
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urel  =  “T  '  “A 


U)j  =  PT  ij  +  Qt  jT  +  R-J-  kT 


“A  "  PA  ^  +  QA  JA  +  RA  kA 


Finally,  the  relative  orientation  of  the  two  coordinate  systems  is 
required  to  allow  coordinate  transformation  of  vectors.  This  may  be  ac¬ 
complished  in  terms  of  the  usual  Euler  angles  0,  <p  ,  and  <p  for  each  vehicle 
where 

<|i  =  inertial  heading  angle 
0  =  angle  of  elevation  above  local  horizontal 
<j>  =  roll  angle 


Let 


T  (0) 


T  (*) 


T  U) 


coso  0  sin© 

0  1  0 

-si  no  0  cos© 

cosi p  -sinty  0 

si  n»p  cosi|»  0 

_0  0  1_ 

“1  0  0 

0  cos<f>  -sin<f> 

0  sin<ji  cos 


-1  T 

Then,  noting  that  T  (a)  =  T  (a)  =  T(-a),  we  can  express  V 
terms  of  their  components  in  the  attacker's  coordinates,  since 


and  wT  in 
T  1 


rr 

*A 

^T 

=  T(-*t)  T(-0t)  T(-*t)  tua)  T(0a)  T(*a) 

Ja 

Ft. 

kA 

.  < 
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After  determining  these  components  we  can  obtain  the  vectors 


V  |  "  UTt  ]A  +  VTt  jA  +  WT  kA 

[  j  I  i  I 


,JT  PT  ’A  +  °T  h  +  RT  KA 


So  we  may  write 


Jrel  "  (PT  '  PA^  ]A  +  (°T  "  V  jA  +  ^RT  ~  V  P'A 


or  u  ,  =  P  ,  i .  +  Q  ,  j .  +  R  ,  Kft 

rel  rel  A  wrel  A  rel  A 


Therefore,  we  obtain  the  resulting  desired  relations  governing  the 
target's  motion  in  attacker's  coordinates  from  Equation  Cl 

UT  =  UT  +  w,  Q  ,  -  VT  R  . 

T  Tj  T  rel  T  rel 

v,  --  vT  +  uT  r  .  -  w,  p  , 

T  Ty  T  rel  T  rel 

WT  =  W,  +  VT  P  ,  -  Ux  Q  , 

T  T-j-  T  rel  T  rel 

Defining  the  above  quantities  in  terms  of  steady  state  and  perturbation 
quantities  yields  the  steady-state  equations 

UT  =  UT  +  WT  Qrel  *  VT  Rrel 

'l  T]  'l  reil  ‘l  re'l 

VT  =  VT  +  UT  R  ,  -  W,  P  , 

T1  TT1  T1  rell  T1  rell 

I 

WT  =  WT  +  VT  P  ,  -  u,  o  , 

T1  TT1  T1  re'l  T1  rell 


Here  note  that  UT  r  VT  =  W  =  0  since  the  steady-state  target  velocity 

T  t  T 

1  1  1 

is  constant  in  the  target  coordinate  system.  The  above  relations  also  show 


the  required  conditions  for  the  target  steady-state  velocity  to  be  constant 


in  the  attacker’s  coordinate  system  as  well. 


These  conditions  are 


^relj  =  VT-j  Rrel  ] 


UT] Rrel 1  =  WT] Pre1 1 


VT1  Prel^  =  ^T^rel  ^ 


Finally,  the  perturbation  equations  for  the  target's  velocity  components 
in  the  attacker's  coordinates  are 


UT  =  uTt  +  \  qrel  +  WT  ^relj 


‘  VT1  rrel  '  VT  Rrel , 


VT  =  vTt  +  UT1  rrel  +  UT  Rrel1 

"WT-jPrel  "  WT  Prel., 

WT  =  wTt  +  VT]  prel  +  VT  Prel1 

"UT1  qrel  "  UT  Qrel., 

Now  if  the  target  is  simply  in  a  steady  flight  condition,  note  that  the  tar- 

+  A 

get's  perturbation  acceleration  terms  uT  =  vT  =  wT  =  0. 

't  't  't 
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APPENDIX  D 

FI 06  PERTURBATION  EQUATIONS 

In  this  section,  the  linearized  equations  for  the  vehicle  dynamics 

★ 

are  summarized.  From  the  reference  the  perturbation  equations  for  a 
general  steady-state  flight  condition  are  as  follows: 


u  -  V-j r  -  R^v  +  W^q  +  Q-jW  =  -gecosG^  +  Xyu  +  Xww  +  X6  6^ 
v  +  Ujr  +  R^u  -  W^p  -  P^w  =  -gesin^  sino^ 

+  g4>cos$i  cosg^  +  Yvv  +  YpP 

+  V  +Ya  5a+  Ya  6p 
r  6A  A  sR  R 

w  -  U^q  -  +  V^p  +  P^v  =  -gecos4>.|  sino^ 

-  g<j>sin«t>^  cosG-j  +  Zuu 

+  Z  w  +  Z  <SC 
w  fig  E 

q  +  Tyy  <‘xx  -  'zz'<Plr  +  RlP)  +  'xz\y  -  2V> 


M  u  +  M  w  +  M*w  +  M  q  +  M.  &r 
u  w  w  q^  6r  E 


P  -  I 


xz/,  '  -  ‘xz/,  (Plq  +  Qlp> 

*xx  *xx 

+  1  <‘zz  •  'yy)(Rlq  *  Qlr)  ■  LS6  *  Lpp 

X  X 


*  Lrr  *  LJA  SA  +  L6r  5R 


f  '  'xz/,  P  +  <V  +  Rlq>  *  Tzz  <‘yy  -  Ixx»Plq  *  V> 


zz 


■  V  *  V  *  V  *  \  SA  +  \  SR 


where  in  these  relations  the  (*)j  quantities  (e.g.,  Uj,  Vj,  W^)  are  constant, 

steady-state  values,  and  X  ,  Z.  etc.  are  the  dimensional  stability  deriva- 

w  w 

* 

Roskam,  J.,  Flight  Dynamics  of  Rigid  and  Elastic  Airplanes,  Lawrence, 

Kansas,  1979. 
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1 


tlves  (see  the  above  reference). 

The  case  to  be  evaluated  consists  of  a  level,  4-g  (n  =  4)  turning 
condition  at  10,000  ft.  altitude,  Mach  *  0.72.  Therefore,  for  this  case, 
the  steady-state  parameters  are  given  as  follows  (equations  developed  In 
the  aircraft  stability  axis) 

e1  =  0. 

U-j  =  V  =  775  ft/sec  cos  ^  =  ^ 

V1  =  W1  =  0.  =  75.52  deg. 

Turn  rate,  $  =  g(tan  *1)/U1  *  .161  rad/sec. 

Roll  rate,  =  0 

Pitch  rate,  Qj  =  f  sin  =  .156  rad/sec 
Yaw  rate,  R1  =  i  cos  =  .040  rad/ sec. 

The  angle  of  attack  required  for  trim  (Ref.  )  is  14.5°  and  the  elevator 
deflection  for  trim  is  -20  deg  at  this  altitude  and  velocity.  The  Inertias 
(in  the  body  axis)  are 

Ixx  =  2.126  x  104  slug- ft2 

5  2 

I  =  2.035  x  103  slug- ft 

I  -  2.175  x  105  slug- ft2 

IX2  *  7.316  x  103  slug-ft2 

The  reference  wing  area  Is  695  ft  and  the  weight  Is  33,000  pounds. 

In  the  stability  axis,  the  Inertias  are 

Ixx  -  2.999  x  104  slug-ft2 

I  »  2.088  x  105  slug-ft2 

zzs 

I  -4.111  x  104  slug-ft2 
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Finally,  the  dimensional  derivatives  are  estimated  to  be  (for  c=m.a.c. 


=  23.8  ft) 

X  =  -.0924  (sec-1) 
u 

X  =  -176.0  (ft/sec2) 
a 

X.  =  -24.4  (ft/sec2) 

6E 

Zu  =  -.437  (sec'1) 

Z5  =  -987.  (ft/sec2) 

Z  =  -288. (ft/sec2) 

6E 

Mu  =  -.00223  (ft-1 sec"1) 

M  =  -11.5  (sec-2) 
a 

M*  =  -.785(. 3485) (sec-1) 
a 

Mq  *  -.785(.3485)(sec-1 ) 

M  =  -17.36  (sec"2) 

6E 

N$  =  6.77  (sec-2) 

Np  =  -.176  (sec"1) 


Ye  =  -177.  (ft/sec2) 
Yp  =  1.76  (ft/sec) 

Yr  =  8.13  (ft/sec) 

Y.  =  84.9  (ft/sec2) 
A 

Y.  =  36.5  (ft/sec2) 
6r 

Lg  =  -47.6  (sec"2) 

Lp  =  -.0613  (sec-1) 

Lr  =  2.77  (sec"1) 

l  =  -57.4  (sec"2) 
6A 

lx  =  .567  (sec-2) 

6R 

Nr  =  -.842  (sec"1) 

N.  =  -5.93  (sec"2) 
A 

N.  =  -3.98  (sec"2) 


Finally,  in  the  steady-state  condition  being  considered,  the  aileron/ 
rudder  interconnect  results  in  the  following  relation 

fiR  "  fiR  commanded  1-1 6A 

In  addition,  the  feel  system  for  this  vehicle  is  estimated  to  provide 
the  following  stick  gains 

6_  =  .0467  6r  (In.)  6  =  .0273  6.  (in.) 

E  Estick  M  Stick 

«R  *  .0738  3^  (1n.) 
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Combining  all  above  parameters,  the  resulting  vehicle  dynamics  in 
this  turning  condition  are  expressed  as 

u  =  .0924  u  -  297. a  -  32.26  +  31.26  -  1.146F 

st 

k  =  -.00036  u  -  1 . 272a  +  q  -  .040$  -  .01736P 

St 

q  =  -.00013  u  -  11.21a  -  .547q  +  .0109* 

+  .0353  p  -  .0162r  -  .8066p 

Stick 

B  «  -.2286  +  .0104*  +  .00227p  -  .99r 

-  . 00005u  +  .001566.  +  . 00347 6„ 

St  Kped 

p  =  - . 0245q  -  77.96  +  . 1 98p  +  5.27r 

-  2.0915.  +  . 6086D 

flst  Kped 

r  =  . 01 27q  +  22.1b  -  .344p  -  1.852r 

+  .3706.  -  ,4136p 

st  Kped 

In  addition  to  the  vehicle  dynamic  model,  the  vehicle  kinematics  are 
described  by  the  perturbation  relations 

psi-y 

q  =  0  cos*j  +  ^  *  cos*i  +  *  sin$j 
r  =  *sin$^  +  *cos#j  -  e  sin$^ 

In  our  case,  solving  for  6,  *,  and  *  yields 

*  =  p  +  .161e 

0  *  -.161*  +  .25q  -  .968r 

*  *  .968q  +  .25r 
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APPENDIX  E 


SYSTEM  STATE  VARIABLE  MODEL 


SYSTEM  DYNAMICS  ARE:  kCOT  =  AX  +  fiU*EW# 
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APPENDIX  F 


MANAGEMENT  SUMMARY 

Summarized  in  this  appendix  is  the  list  of  written  publications, 
including  those  planned;  data  on  the  advanced  degree  student  supported 
by  this  research;  and  the  summary  of  conference  papers  presented  and 
interactions  with  the  Air  Force  Wright  Aeronautical  Laboratories  (formerly 
the  Flight  Dynamics  Lab). 

Publications  include  the  following: 

1)  D.K.  Schmidt,  "Pilot-Optimal  Augmentation  for  the  Air-to-Air 
Tracking  Task,"  to  appear  in  the  AIAA  Journal  of  Guidance  and 
Control ,  1980. 

2)  S.N.  Prasad  and  D.K.  Schmidt,  "Closed-Loop  Pilot/Vehicle  Analysis  for 
Multi-Axis  Air-to-Air  Tracking,"  in  preparation  for  submission  to  the 
AIAA  Journal  of  Guidance  and  Control. 

One  graduate  student  was  supported  for  one  semester  and  one  summer 
session  on  this  program  and  much  of  the  technical  results  herein  are  in¬ 
cluded  in  his  thesis.  This  student,  Mr.  S.N.  Prasad,  will  receive  the 
Master's  Degree  in  Aeronautical  and  Astronautical  Engineering  from  Purdue 
University  in  August,  1980,  and  his  M.S.  thesis  will  be  published  at  that 
time. 

Other  dissemination  of  this  studies  results  has  been  accomplished  via 
papers  presented  at  technical  conferences.  These  include  presentation  of 
the  first  paper  cited  above  at  the  1979  AIAA  Guidance  and  Control  Conference 
held  at  Boulder,  Colorado,  August  6-8,  1979.  In  addition,  a  presentation 
entitled  "Multi-Axis  Tracking  Via  an  Optimal  Control  Pilot  Model,"  was  made 
at  the  16th  Annual  Conference  on  Manual  Control  at  Cambridge,  Mass,  in 
May,  1980. 
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Finally,  throughout  the  course  of  this  program,  frequent  contacts 
have  been  made  with  the  members  of  the  technical  staff  of  the  Flight 
Control  Division,  Air  Force  Wright  Aeronautical  Laboratories  (formerly 
the  Flight  Dynamics  Laboratory).  These  contacts  have  included  partici¬ 
pation  in  workshops  (1979  Flying  Qualities  Symposium,  October,  1979,  and 
the  1980  Stability  and  Control  Workshop,  April,  1980)  held  at  Wright- 
Patterson  AFB,  Ohio.  Also,  informal  telephone  contacts  with  Mr.  Ronald 
0.  Anderson  and  Mr.  Frank  George  of  the  Wright  Aeronautical  Laboratories 
have  been  made  a  number  of  times  throughout  this  study. 


